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EXECUTIVE  SUMMARY 


Contract  F33615-01-D-2103  was  awarded  on  19  April  2001  as  an  Indefinite  Delivery, 
Indefinite  Quantity  (IDIQ)  contract  with  a  ceiling  of  $13,553,119  and  an  anticipated  expiration 
date  of  3 1  October  2004  or  as  extended.  Multiple  delivery  orders  were  expected  to  be  awarded 
under  the  contract.  Delivery  Order  no.  2,  the  subject  of  this  final  report,  was  awarded  on 
September  24,  2001  and  eventually  funded  to  April  30,  2003  with  amendments. 

The  effort  funded  by  this  contract  encompassed  two  principal  objectives: 

1 .  The  continuation  of  a  project  funded  by  Contract  F33615-00-C-2030,  which  created  a  silicon 
carbide  semiconductor  device  prototyping  capability  at  Mississippi  State  University  known 
as  the  Mississippi  Center  for  Advanced  Semiconductor  Prototyping  (MCASP). 

2.  The  extension  of  the  earlier  work,  which  demonstrated  a  power  Schottky  barrier  diode  (SBD) 
prototyped  in  pre-production  lot  size  (i.e.,  many  more  parts  than  is  typical  of  research),  to  a 
unipolar  silicon  carbide  semiconductor  power  switch,  the  Vertical  Junction  Field  Effect 
Transistor  (VJFET). 

3.  The  continued  involvement  by  a  direct  commercialization  partner,  namely  SemiSouth 
Laboratories,  Inc.,  in  the  development  of  a  realistic  path  to  system  insertion  for  the 
technology  funded  by  this  contract  and  its  predecessor  (F33615-00-C-2030),  via  a 
public/private  partnership  involving  technology  licensing  and  subcontracting. 

The  two  larger  objectives,  the  development  of  a  viable  power  semiconductor  switch  and 
the  creation  of  a  realistic  commercial  path  to  DoD  system  integrators  (the  large  defense  prime 
contractors)  cannot  be  accomplished  in  one  delivery  order.  This  Delivery  Order  no.  2  final  report 
examines  another  aspect  of  the  project,  namely,  the  design  of  a  viable  power  vertical  junction 
field  effect  transistor.  A  design  based  on  a  self-aligned,  gate-implanted,  trenched  source-gate 
junction  FET  was  selected  for  its  near  term  technological  readiness  and  it  long  term 
manufacturability.  Delivery  Order  no.  2  concentrated  on  several  key  aspects: 

•  Optimized  design  of  the  fundamental  static  induction  transistor  that  is  the  heart  of  the  power 
VJFET.  This  includes  critical  dimensions  of  the  gate  and  the  channel  fingers. 

•  Optimized  design  of  an  appropriate  edge  termination  suitable  for  scaling  the  gate-drain  diode 
of  the  VJFET,  where  the  drain  potential  is  blocked,  to  high  voltages  in  excess  of  1.2  kV. 

This  final  report  will  focus  on  the  design  of  the  critical  static  induction  transistor.  An 
appendix  has  been  added  that  documents  the  edge  termination  design  optimization  effort  in  the 
form  of  a  paper  published  at  the  International  Conference  on  Silicon  Carbide  and  Related 
Materials  (ICSCRM)  held  in  Lyon,  France  in  2003.  However,  that  paper  documents  the  first 
investigation.  Continuing  investigation  in  later  delivery  orders  has  led  to  further  understanding 
of  the  issues,  including  a  need  to  improve  the  edge  termination  and  passivation  approach  applied 
to  the  power  trench  VJFET.  It  is  now  recognized  that  an  integrated  approach  using  features  of 
the  gate-drain  pn  diode  lead  to  the  best  results,  which  will  be  documented  in  later  reports. 

The  follow-on  delivery  orders  (3-5)  will  report  the  performance  of  a  critical  transistor 
component  (ohmic  contacts  -  Delivery  Order  no.  3)  and  the  performance  of  the  numbered 
VJFET  developmental  lots  (Delivery  Order  nos.  4  and  5).  A  single  final  report  will  be  issued 
covering  both  of  the  later  delivery  orders.  Delivery  order  6  concludes  the  contract,  and  will 
summarize  the  final  results  of  device  development,  applications  development,  and  the  initial 
reliability  testing  studies.  A  single  report  will  be  issued  covering  Delivery  order  6,  including  the 
enhancement. 
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CHAPTER  1 
INTRODUCTION 

1.  Background  :  Review 

Silicon  Carbide  (SiC)  as  a  semiconductor  material  is  provided  with  several  special 
properties.  The  primary  research  areas  for  SiC  have  been  high  power  RF  devices,  high 
temperature  electronics  and  radiation  hard  electronics  (space  electronics  and  nuclear  reactor). 
This  can  be  attributed  to  special  material  properties  of  SiC.  The  large  bandgap  of  SiC  (3.2eV  for 
4H  polytype)  and  high  electron  field  strength  (3.0  MV/cm)  allows  for  high  voltage  applications 
and  reduces  series  resistance  in  the  lightly-doped  drift  region.  The  high  saturated  drift  velocity  of 
SiC  is  also  advantageous  in  the  development  of  high  speed  devices  such  as  high  frequency 
amplifiers  [1],  Table  1-1  shows  a  comparison  of  material  parameters  for  various  semiconductors. 

SiC’s  Static  Induction  Transistor  (SIT)  is  well-suited  for  high  power  and  high  frequency 
(71  kW  demonstrated  at  UHF).  The  leading  parameters  for  design  are  developed  through  careful 
modeling  by  means  of  simulations  and  inclusion  of  all  significant  device  physics,  this  thesis 
demonstrates  that  good  agreement  is  reached  between  theoretical  prediction  and  results  of 
various  simulations  ,  such  as  temperature,  saturation  drift  velocity,  and  contact  resistance. 

The  SiC  Static  Induction  Transistor  is  based  on  the  original  SIT  concepts  developed  in 
1950  [2]  by  Nishizawa,  which  have  received  continued  development  for  the  past  5  decades  [3-5]. 
As  compared  to  a  Si  Bipolar  Junction  Transistor  (BJT),  the  SIT  has  advantages  in  being  a 
majority  carrier  device,  voltage-controlled,  and  often  higher  breakdown  voltage  and  input 
impedance.  While  most  SIT  development  has  occurred  in  Si,  renewed  interest  in  the  use  of  SiC 
SIT’s  over  the  past  10  years  has  occurred  because  of  SiC’s  ten- fold  advantage  (3  MV/cm)  in 
critical  field  strength  and  nearly  two-field  advantage  in  saturated  electron  velocity 
(l.6  ~  2.2xl07cm/s)  [6-10].  Additionally,  SiC  has  an  extremely  high  thermal  conductivity 
(0k  =3  W  /cm- k)  which  far  surpasses  that  of  Si  and  GaAs,  and  allows  an  increased  junction 
operating  temperature  (7)  )to  nearly  400°C  because  of  the  high  0k  and  dramatically  reduced  pn 
junction  leakage  currents  resulting  from  its  wide  (E  =  3.2  eV  for  the  4H-SiC  polytype) 

bandgap.  By  use  of  SiC,  it  is  thus  possible  to  increase  the  voltage  (power)  rating  of  high- 
frequency  amplifiers,  and  also  increase  the  power  density  by  a  factor  of  4  to  10  in  UHF  to  S- 
band  transistors.  Parasitic  capacitances  are  reduced  by  shrinking  the  device  area  while 
maintaining  the  power  level  and  higher  input/output  impedance  of  the  devices  often  requires 
little  or  no  impedance  matching.  Limitations,  which  have  slowed  commercial  insertion  of  SiC 
devices  to  date,  include  wafer  size  (75  mm  maximum  diameter  available  in  2002),  wafer  cost, 
device/epitaxy  processing,  and  wafer  defects.  Currently,  solid  state  devices  are  available  provide 
a  total  power  of  tens  of  watts  at  1  GHz,  such  as  the  silicon  based  bipolar  transistor,  and  a  few 
watts  at  10  GHz,  like  GaAs  based  MESFET.  When  higher  power  is  needed,  one  needs  to  either 
combine  the  output  power  of  multiple  devices  or  use  vacuum  tubes,  which  are  still  uncontested 
at  very  high  power  levels,  capable  of  up  to  a  few  hundred  kilowatts  at  5  GHz  [11],  High  total 
power  for  a  single  device  means  that  a  smaller  number  of  devices  are  needed  to  build  a  high 
power  amplifier.  This  should  lead  to  higher  reliability,  easier  circuit  design  and  lower  cost. 
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2.  Properties  of  SiC 
2.1  Crystal  Structure 

SiC  exhibit  a  one-dimensional  polymorphism  called  polytypism.  An  almost  infinite  number 
of  SiC  polytypes  are  possible,  and  approximately  200  polytypes  have  already  been  discovered. 

[12]  SiC  polytypes  differentiated  by  the  stacking  sequence  of  each  tetrahedrally  bonded  Si-C 
bilayer.  With  the  exception  of  2H  and  3C,  all  of  the  polytypes  form  1-D  superlattice  structure 

[13] .  The  polytypes  are  divided  into  three  basic  crystallographic  categories;  cubic  (C),  hexagonal 
(H),  and  rhombohedral  (R).  Cubic  SiC  has  only  one  possible  polytype,  and  is  referred  to  as  3C- 
SiC  or  P-SiC.  Each  SiC  bilayer  can  be  oriented  into  only  three  possible  positions  with  respect  to 
the  lattice  while  the  tetrahedral  bonding  is  maintained.  If  these  three  layers  are  arbitrarily 
denoted  A,  B,  and  C,  and  the  stacking  sequence  is  ABC/ABC...  then  the  crystallographic 
structure  is  cubic  zinc  blende  and  referred  to  as  3C-SiC.  It  possesses  the  smallest  bandgap 
(~2.4eV)  [14],  and  one  of  the  largest  electron  mobility  of  all  the  SiC  polytypes,  and  has  been 
grown  on  6H-SiC  substrate.  If  the  stacking  of  the  bilayers  is  ABC/ABC...,  then  referred  to  as  2H- 
SiC,  stacking  sequence  with  ABC/ABCB... and  ABCACB/ABCACB ...,  is  referred  to  as  4H-SiC 
and  6H-SiC.  4H-SiC  consists  of  equal  number  of  cubic  and  hexagonal  bonds.  6H-SiC  is 
composed  of  two-thirds  cubic  bonds  and  one-third  hexagonal  bonds. 

These  SiC  polytypes  have  substantially  different  thermal,  optical,  and  electrical  properties 
due  to  their  different  crystal  structure.  The  properties  of  the  commonly  used  SiC  polytypes, 
which  are  most  important  for  device  applications,  are  comparing  with  the  similar  properties  of 
silicon  and  gallium  arsenide  in  Table  1.1. 

The  properties  of  SiC,  such  as  wide  band  gap,  excellent  thermal  conductivity,  high 
breakdown  field,  high  saturation  drift  velocity,  provide  a  great  deal  of  attractiveness  in  high 
speed,  high  power  devices  working  at  elevated  temperatures  and  radiation,  such  as  high  power 
switch  and  high  power  amplifier.  Also,  SiC’s  properties  of  high  chemical  stability  and  radiation 
tolerance  have  attraction  for  aviation,  nuclear,  and  space  applications. 


Table  1.1  Properties  of  SiC  polytypes 


Property 

Si 

GaAs 

4H-SiC 

6H-SiC 

. 

3C-SiC 

Band  Gap(eV) 

1.12 

1.42 

3.2 

3 

3.2 

Breakdown  Field 
(MV/cm) 

0.6 

0.6 

3.0 

3.0 

3.0 

Saturated  Drift 

Velocity  (cm/s) 

107 

107 

2.2xl07* 

1.9xl07* 

2.5xl07 

Thermal  Conductivity 
(w/cm  °C) 

1.5 

0.5 

4.9 

4.9 

5 

Electron  Mobility 
(cmVVs)* 

1100 

6000 

900 

370 

750 

Hole  Mobility 
(cm2/Vs)+ 

420 

320 

115 

90 

40 

at  doped  concentration  of  10 17  cm'3 
*:[15] 
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2.2  Energy  Band  Structure 


The  energy  band  structure  of  SiC  strongly  depends  on  the  polytypes.  Both  experimental  and 
theoretical  results  show  that  the  valence  band  maxima  for  all  common  polytypes  are  located  at 
the  center  of  the  Brillouin  zone  ( T  point).  For  all  polytypes  but  2H-SiC  the  conduction  band 
minimum  of  4H-SiC  is  found  to  be  at  the  M  point  of  the  Brillouin  zone.  The  indirect  gap  varies 
with  the  polytype  reaching  values  between  Eg  =  2.3  eV for  3C-SiC  and  Eg  =  3.2  eV for  4H-SiC 

has  the  largest  energy  band  gap  among  all  common  SiC  polytypes. 

The  temperature  dependent  energy  band  gap  is  not  known  for  4H-SiC.  Experiment 
results  Ref.  in  [16]  show  that  the  temperature  coefficients  of  energy  band  gap  for  6H-SiC  and 
3C-SiC  are  -  3.3  x  10  4 eV/K  and  - 5.8  x  10  4  eV/K,  respectively,  Since  crystal  structure  of  4H- 
SiC  is  closer  to  6H-SiC  than  to  3C-SiC,  following  the  same  approach  as  in  Ref.  [17]  and  using 
the  temperature  coefficient  of  the  energy  band  gap  of  4H-SiC,  the  following  temperature 
dependent  band  gap  of  4H-SiC  is  obtained 

Eg  =3.19-3.3xl0~4(f-300/f)  eV  (1.1) 

2.3  Intrinsic  carrier  concentration 

In  semiconductor  4H-SiC,  effective  densities  of  states  in  the  conduction  and  valence  bands  can 
be  expressed  as 

4  7mckt 
h2 

47m  vkt 
h 2 


(1.2) 

(1.3) 


where  k  is  the  Boltzman  constant.  For  4H-SiC  mc  (300/i)  =  0.16m0  mv  (300£)«1.2m0  in  room 
temperature.  The  intrinsic  carrier  concentration  at  thermal  equilibrium  is  given  by 
, -  f-Ee  (rVi 

/  A  T  A  T  _  £  ^  '  \  /  1  A  \ 


n,  =  JNCNV  exp 


The  temperature  dependent  intrinsic  carrier  concentrations  of  4H-SiC,  Si,  and  GaAs  are  shown 
in  Figure  1.1. 
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Figure  1.1  Intrinsic  carrier  concentrations  vs  temperature  of  4H-SiC,  Si  and  GaAs. 

(calculated  from  Eq.  1.4) 


2.4  Low  field  mobility 

The  low  field  mobility  model  for  4H-SiC  can  be  expressed  as 

delta  /  fj. j  \czn, 

ju  =  jumm  +- 

r^n,p  r~n,p 


1  + 


f  \/n,p 

ND+N/ 


V 


Nu 

n,p  J 


300 K 


(1.5) 


where  the  electron  set  of  data  are  ju™'"  =  0  cm 2  /(V  sec),  judnelta  =  947  cm 2  /(V  sec), 
=  1.94x  1017  cm2  /(F  sec),  yn  =  0.61  and  an  =  -2.15 ;  the  hole  set  of  data  are 


//“n  =  15.9  cm2  /(V sec),  judnelta  =  108.1  cm2  /(V sec),  N%  =1.76xl019  cm2  /(v sec), 
yp  =  0.34  and  ap  =  -2.15 . 

Note:  this  mobility  model  does  not  contain  the  effects  of  electron-hole  scattering,  which  may 
be  important  in  power  devices  and  operating  under  high- injection  conditions  [17]. 
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Figure  1.2  Low  field  electron  and  hole  mobility  of  4H-SiC,  Si,  and  GaAs  at 
room  temperature 

The  low  field  mobility  of  4H-SiC  higher  than  those  of  6H-SiC  for  both  electrons  and  holes; 
the  mobility  of  Si  is  higher  than  that  of  4H-SiC,  and  that  of  GaAs  is  highest.  Thus  GaAs  has 
tremendous  advantage  in  high  speed  application  field  over  both  Si  and  4H-SiC  at  low  fields.  The 
low  field  mobility  of  4H-SiC,  Si,  and  GaAs  at  room  temperature  increases,  the  lattice  atoms 
vibrate  more  about  their  mean  position,  thus  effectively  increasing  in  “size”  and  thereby 
increasing  the  chance  of  collisions.  The  mobility  therefore  decreases  with  increase  T.  The 

-3 

theoretical  dependence  is  //go T  2  [18]  and  it’s  shown  in  Figure  1.3. 

2.5  Thermal  conductivity 

The  other  important  physical  property  of  SiC  is  its  high  thermal  conductivity.  Typically,  the 
thermal  conductivity  depends  on  polytypes  and  doping  but  the  thermal  conductivity  of  SiC 
exceeds  that  of  copper  (4.01  Wcm~lK~l),  silver  i^A29Wcm~x ,  Al203,  and  is  about  5  times 

higher  than  that  of  Si.  This  is  mainly  due  to  the  phonon  dispersion  relation  distribution.  Also, 
this  property  leads  to  a  single  design  of  heat  dissipation  of  the  device  and  the  circuit  even  in  high 
temperature  and  high  power  [19]. 
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Figure  1.3  Mobility  vs  doping  and  temperature  for  electron  and  holes  4H-SiC, 
GaAS,  Si. 
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CHAPTER  II 

STATIC  INDUCTION  TRANSISTOR 

1.  Device  Structure 

Both  recessed  gate  and  planar  gate  SIT's  have  been  fabricated  in  Si  and  SiC.  However 
because  of  difficulty  in  deep  ion  implantation  (p-type  dopant  A1  implanted  at  energy  of  380  keV 
has  a  peak  range  of  only  0.42p«t  into  SiC)  and  the  lack  of  diffusion  at  temperatures  below  1800 
°C,  the  recessed  gate  structure  has  been  favored  in  SiC.  Typical  recessed  gate  SIT  structures  are 
shown  in  Fig.  2.1,  with  primary  differences  found  in  the  gates.  The  gate  may  be  a  Recessed 
Gate  -  Bottom  contact  (RG-B),  as  shown  in  Fig.  2.1,  or  a  RG  with  Sidewall  and  Bottom  contact 
(RG-SF3),  as  shown  in  Fig.  2.2  which  provides  a  longer  gate  length  (Lg).  Increasing  Lg 

improves  gate  control,  voltage  gain  (W.  and  transconductance  ( gm ),  while  reducing  the  output 

resistance  (rDS).  The  gate  material  may  be  formed  from  either  a  Schottky  metal  or  a  pn  junction, 
and  can  be  thought  of  as  analogous  to  a  vertical  dual-gate  MESFET  or  JFET  structure, 
respectively.  The  gate  depletes  the  channel  from  each  side  a  distance  (xa)  into  the  channel.  The 
total  channel  thickness  (2a),  trench  width  ((),  and  drift  region  length  (Lgd)  are  illustrated  in  Fig. 
2.1.  Planar  gates  can  also  be  fashioned  in  SiC  [1],  but  the  recessed  gate  with  a  Schottky  barrier 
offers  the  lowest  parasitics  and  most  efficient  gating  [2],  A  PN  junction  gate  can  be  used  with 
higher  gate  resistance  (Rg),  but  offers  less  gate-to-source  leakage  at  a  higher  maximum  junction 
temperature  (T/max )  and  potentially  higher  power  density.  The  PN  gate  leakage  at  high 

temperatures  in  SiC  devices  is  substantially  reduced  (proportional  to  the  intrinsic  carrier 
concentration  of  ni  or  ni2 )  compared  to  Si  or  GaAs  devices,  and  degradation  of  the  Schottky 
contact  is  not  limiting  reliability  if  one  uses  a  PN  gate,  as  shown  in  Fig.2.3.  PN  gate  adds  some 
cost  to  processing  since  additional  p-type  implant  and  activation  anneal  are  required. 
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Figure  2.1  Structure  of  recessed  gate  static  induction  transistor. 
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Figure  2.2  Structure  of  well-suited  recessed  gate  static  induction  transistor 


Figure  2.3  Structure  of  surface  gate  static  induction  transistor. 


2.  Operation 

With  Vg  =  VS  =0  V,  the  conventional  SIT  has  its  channel  completely  depleted  (xd  >  a). 
If  xd  >  a,  which  often  is  done  to  increase  maximum  total  current  and  power,  then  the  SIT  is 

referred  to  as  mixed-mode  [4],  Other  terms  for  the  SIT  include  gridstor,  VFET  (Vertical  FET), 
or  Permeable  Base  Transistor  (PBT)  [15].  PBT’s  often  have  extremely  short  gates 
( Lg  <  a  <  xd ) .  The  VFET  nomenclature  is  generally  reserved  for  devices  with  saturated  drain- 

source  output  curves,  unlike  the  triode  characteristics  of  the  SIT.  The  SIT  name  itself  refers  to 
the  fact  that  the  output  characteristics  are  modified  by  static  induction  and  not  just  the  input 
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(gate)  voltage,  as  is  the  case  for  conventional  FET’s  in  the  saturation  mode.  There  are  three  basic 
conduction  mechanisms  possible  when  the  SIT  is  on,  known  as  the  ohmic,  exponential,  and 
Space  Charge  Limited  Current  (SCLC)  regions,  which  are  identified  in  Fig.  2.4  for  a  typical  RG- 
SB  SIT  in  SiC  with  a  channel  and  drift  region  doped  at  5xl015cm“3.  (One  may  increase  the 
channel  doping  to  increase  the  total  maximum  current.)  The  transistor  conduction  mechanisms 
are  two-dimensional,  and  a  short  basic  description  of  the  various  regimes  is  given.  Current  flow 
is  modified  by  gate,  drain,  and  source  potentials,  through  a  saddle  point  minima  between  the  gate 
electrodes,  as  shown  in  Fig.  2.5.  where  the  minima  is  located  at  the  mid-point  of  the  channel. 
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Figure  2,4  SIT’s  three  regions  of  current  flow. 
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Figure  2.5  Conduction  Band  Energy  of  one  half  drift  channel. 
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Figure  2.6  Conduction  Band  Energy  in  some  bias. 
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Figure  2.7  Conduction  Band  Energy  in  the  center  of  channel  (X=0)  for 
different  drain  voltages. 


2.1  Ohmic  region 

In  the  ohmic  (linear)  region,  the  channel  is  not  completely  pinched  off  by  the  depletion 
regions  under  fairly  low  gate  and  drain  bias  (see  Fig.  2.5,  Fig.  2.6,  and  Fig.  2.7).  The  current 
from  drain-to-source  (I ds)  is  controlled  by  ohmic  resistances  found  between  the  depletion 

regions  under  (beside)  the  gate,  undepleted  portions  of  the  drift  region,  and  the  partially 
depleted  portion  of  the  channel  immediately  above  the  gate.  The  limiting  portion  is  of  course 
the  region  between  gate  depletion  regions,  which  is  often  very  narrow.  Each  channel  is 
depleted  from  the  gate  on  both  sides  by  xj  which  is  given  by  Eqn.  (2.1),  where  e  s=  10.0  for 
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SiC,  s  0=  8.854  x  10'14  F/cm 2  and  No  is  the  channel  doping.  Vbi  and  V  represent  the  built-in 


potential  and  gate  voltage,  respectively,  s 

_  |2Vo( Vu-V,) 


where:  Vbi  =</>B  — 

q 

The  quantity  (j)B  is  the  Schottky  barrier  height  (typically  1.1  eFfor  Ni  on  SiC)  at  the  metal- 
SiC  junction.  The  current  density  in  this  region  can  be  given  by  JD  =  q-vnD  • ND ,  and  the 
area  is  simply  A=2.A  eff  .  Z  per  finger.  The  quantities  Aej^and  Z  are  the  undepleted  channel 

width  and  channel  depth,  respectively.  In  this  low-field  region,  the  current  is  not  limited  by 
vsat,  but  rather  proportional  to  the  effective  open  area,  channel  doping,  and  low-field 

electron  velocity  ( vnD ).  For  this  reason,  the  4H  polytype  of  SiC  is  the  commercially- 
available  preferred  polytype  because  it  has  the  highest  electron  mobility. 
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2.2  Exponential  region 

In  the  exponential  region,  even  though  the  channel  is  depleted  with  a  low  gate  bias,  a 
high  drain  bias  influences  conduction  over  the  potential  barrier  in  the  channel  from  drain  to 
source.  The  potential  barrier  in  the  channel  (as  shown  in  Fig.  2.5  and  2.6  for  half  of  a  1.2 
\xm  channel  thickness)  between  the  two  gate  electrodes  has  barrier  lowering  in  the  mid-point 
region  exponentially  dependent  upon  the  high  drain  bias.  Since  the  exponential  relationship 
exists,  the  current  in  this  region  can  be  described  below  in  Eqn.  (2.2)  as  exponentially 
dependent  upon  the  amount  of  barrier  lowering  in  the  channel  (AO),  which  is  a  function  of 
the  drain  and  gate  bias  [6]. 

/  =  V"r  (2-2) 

where  q  is  the  electron  charge,  k  is  Boltzman’s  constant,  T  is  the  absolute  temperature,  and 
1 0  is  a  constant  that  depends  on  the  channel  geometry. 


2.3  Space  Charge  Limited  Conduction  region 

When  the  current  density  becomes  high,  a  large  drain  voltage  induces  a  large 
number  of  carriers  into  the  gated  and  drift  region.  Also,  when  the  carrier  concentration 
becomes  higher  than  the  doping,  the  current  becomes  space  charge  limited.  One  can  modify 
the  treatment  by  Mott  and  Gurney  [7,8]  and  arrive  at 


,  V 

J D  —  2.£s£0VSAt  ^2 


(2.3) 


where,  V  is  the  voltage  drop  from  the  contact  (drain  terminal)  to  the  saddle  point  minima  and 
L  is  the  distance  from  the  drain  to  the  saddle  point  minima  (~  Lsg  +  Lg  +  Lgd).  v  SAT  is  the 


saturated  electron  drift  velocity.  Eqn.  (2.3)  can  be  expressed  in  terms  of  total  current  [2]  as 
shown  in  Eqn.  (2.4).  Note  that  total  current  for  similar  voltage  and  frequency  rated  Si  and 
SiC  devices  will  be  quite  different  as  Lgd  would  be  an  order  of  magnitude  thinner  and  v  SAT 


would  be  1.5  to  2X  higher  in  the  case  of  SiC.  However,  in  using  these  equations,  there  is 
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an  intrinsic  assumption  made  regarding  the  SIT  operation,  which  is  explicitly  that  the 
electric  field  remains  high  enough  (DlxlO5  V/cm  in  SiC)  throughout  the  device,  so  that  the 
electron  velocity  remains  in  saturation.  While  this  is  a  simplifying  assumption,  which  makes 
analysis  easier,  it  can  easily  be  violated,  and  can  have  significant  impact  on  device  operation 
(see  discussion  regarding  temperature  limits  below  for  example). 


1(1  ~  ^£s£0Vsat 


V ■ 


,2  ^ eff 


(2.4) 


3.  Limitations  for  Design 

3.1  Thermal  Limit  on  Total  Power 

Most  power  devices  have  a  heat  sink  that  is  usually  attached  to  the  backside  of  the 
wafer.  The  generated  heat  must  be  transferred  to  it.  Assuming  a  45°  spreading  angle,  the 
associated  thermal  resistance  can  be  computed  by 


R 


TH 


,Jtm+W)z) 

(T„(  +  Z)W) 


(2.5) 


where,  W  and  Z  are  width  and  depth  of  device,  and,  for  a  given  power  density  per  area  P” 
we  find  the  rise  in  device  temperature  AT  for  steady  state  conditions. 


AT  =  RthP" 


(2.6) 


The  thermal  resistance  Rm  has  a  sizable  spreading  part  that  scales  sublinearly 
with  an  increase  of  the  device  area,  A-WZ  This  means  an  increase  in  internal  device 
temperature  for  larger  area  devices  and  a  give  temperature  of  the  heat  sink.  Typically, 
the  pitch  between  adjacent  mesas,  t+2a,  affects  the  thermal  resistance  as  well  as  the 
thermal  power  density  generated  on  the  top  side  of  the  wafer.  A  larger  pitch  makes  it 
easier  to  dissipate  the  heat  from  the  device.  Additionally  possibilities  to  improve  the 
heat  dissipation  are  a  reduction  of  the  generated  thermal  power  by  a  high  efficiency, 
and  reducing  the  thermal  resistance  by  thinning  the  substrate.  However,  a  thinner 
substrate  may  in  practice  not  improve  the  heat  dissipation  much  because  the  thermal 
conductivity  of  the  SiC  substrate  is  comparable  to  the  thermal  conductivity  of  the 
metals  that  one  would  use  for  a  heat  sink  attached  to  the  backside  of  the  wafer. 

3.2  Impedance  and  Area  Limit  on  Total  Power 

The  output  and  input  impedances  can  use  in  a  microwave  system  have  a  lower 
limit  that  is  fixed  by  circuit  considerations  and  rarely  be  influence  by  the  device  designer. 
The  impedance  at  the  input  and  output  of  a  device  decreases  with  increasing  device 
periphery,  thus  impedance  imposes  a  limit  on  the  total  device  size.  The  input  impedance  of 
a  typical  SIT  is  mainly  capacitive.  For  example,  the  input  capacitance  with  the  output 
shorted  is  Csin,  the  input  impedance  with  output  connected  to  load  is  Cfn.  For  all  cases  of 

interest  Cfn  is  larger  than  Cfn  due  to  the  Miller  effect  [9],  In  Power  microwave,  the  loaded 

input  capacitance;  C  fn  must  be  small  enough  so  that  the  desired  bandwidth  can  be  achieved. 

The  input  impedance  of  the  device  is  realized  by  synthesizing  a  transmission  line  such  that 
the  input  capacitance  of  each  transistor  stage  is  absorbed  into  the  line  impedance.  This 
transmission  line  will  have  a  cut-off  frequency  given  by  [10] 
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CO 


2 


(2.7) 


c 


and, 


Eqn.  (2. 7)  and  Eqn.  (2.8)  lead  to  the  frequency  limit  for  the  amplifier  given  by 

/ 

1 


fc 


,nZ;Ci 


(2.8) 


(2.9) 


The  load  line  is  given  in  Fig.  2.4  between  VBR(lds  =  0)  to  7max (Vds  =  Vknee),  where  VSR  is  the 
static  breakdown  voltage  between  drain  and  source.  7max  is  the  maximum  total  current  in  the 
ohmic  region  at  Vg  =  0  before  conduction  begins  to  be  affected  by  SCLC  and  series 
resistance  denoted  at  Vds  =  Vknee .  The  very  high  VBr  of  SiC  SIT's  (up  to  450  V  for  S-band  parts 
has  been  demonstrated)  allows  the  theoretical  power  density  (PD)  described  in  Eqn.  (2.10)  to 
be  much  higher  than  in  Si.  Since  PD  is  proportional  to  VBR2,  and  SiC's  FeRis  nominally  10X 
that  of  Si  for  similar-rated  parts,  this  translates  to  a  100X  higher  theoretical  power  density 
advantage  for  SiC  SIT's  over  Si.  Achieving  those  advantages  will  also  require  exceptional 
packaging  technology  to  be  implemented  for  SiC  parts. 

1  /  \  (V  -V  Y 

P  =  —  /  (V  -V  1  =  V  BR _ kneej_  f  2  1  Of 

1  D  g  1  max  v  BR  y  knee  J  g^ 

Using  SiC  provides  another  advantage  in  broadband  amplifiers,  which  require  minimum 
output  impedance  from  the  transistor.  Since  scaling  the  area  up  in  a  device  decreases  both 
input  and  output  impedance,  the  area  of  the  device  can  be  constrained  by  this  requirement. 
SiC’s  ability  to  operate  at  higher  voltages  can  increase  the  total  output  power  while  still 
maintaining  some  required  output  impedance. 
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CHAPTER  III 

STATIC  CHRACTERISTICS 


1.  Device  Physics 

The  static  I-V  characteristics  of  the  basic  SIT  were  computed  with  the  three  partial 
differential  equations  (Eqn.  (3-1),  (3-2),  and  (3-3))  self-consistently  for  the  electrostatic 
potential  T  and  for  the  electron  concentrations  n,  respectively.  In  the  SIT,  the  electrical 
behavior  is  governed  by  Poisson's  equation. 

V2vP  =-lLn  +  N+D)  (3-1) 

£ 

and  Continuity  equation  for  electrons. 

V-J  =  Un  (3-2) 

where  \|/,  s,  and  n  are  the  electrostatic  potential,  dielectric  constant,  and  electron  density, 
respectively.  N+D ,  J ,  Un  are  the  concentration  of  ionized  donors,  current  density,  and  net  rate  of 

electron  recombination.  Since  the  SIT  is  a  majority  carrier  device  and  SiC  is  large  bandgap 
material,  the  effect  of  holes  has  been  omitted  in  this  analysis.  The  whole  concentration  will  be 
smaller  that  the  numerical  limit  of  the  simulation  tools.  Carrier  transport  is  given  by  the  drift 
diffusion  equation 

J  =qM„Enn  +  qDnVn  (3-3) 

where  //„ ,  En ,  and  Dn  are  the  electron  mobility,  the  electric  field,  and  electron  diffusivities, 
respectively. 

2.  Voltage  Gain  and  Maximum  Current  Physics 

The  voltage  gain  p  is  a  measure  for  the  possible  large  signal  gain.  It  changes 
considerably  with  drain  voltage,  but  for  the  scope  of  this  work  p  is  defined  for  a  drain  voltage  of 
100  V.  For  practical  devices  the  voltage  gain  should  be  at  around  10,  which  would  correspond  to 
a  device  that  required  -1  OF  on  the  gate  terminal  to  keep  the  device  off  at  100F  on  the  drain.  A 
graphical  interpretation  of  these  parameters  is  shown  in  Fig.  3.1. 


18 


Figure  3.1  On  resistance,  voltage  gain,  maximum  current  and  knee 
voltage. 


The  zero  bias  on-resistance,  Ron  is  simply  the  initial  slope  of  the  drain  current  versus 
drain  voltage  curve  for  zero  gate  voltage 
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<sv^ 
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vd= 0 
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(3-4) 


In  Eqn.  (3-4),  a  small  Ron  is  desirable  for  higher  total  current  and  smaller  knee  voltage. 
The  definition  illustrated  in  Fig.  3.1  places  the  point  of  maximum  current  where  the 


straight  line  through  the  origin  given  by  7max 


intersects  with  the  drain  current  curve 


for  Vg  =  0 .  Since  the  SIT  shows  non-saturating  I-V  curves,  /max  is  not  as  obvious  as  for 
MESFET  type  device.  Also,  should  mention  that  Vgcan  be  up  to  +2.5  V  in  wide  bandgap 

(3.2  eV)  SiC  pn-gate  SIT’s  without  getting  forward  bias  injection  from  the  gate.  But  higher 
7 max  wiH  allow  for  a  higher  power  density  per  unit  width  for  a  given  breakdown  voltage. 
The  knee  voltage  Vknee  is  the  drain  voltage  that  corresponds  to  the  point  of  maximum 

current.  The  knee  voltage  is  important  because  it  is  the  main  factor  that  degrades  the  drain 
efficiency,  which  is  given  by  [1] 


Vd  = 


0.5 


1  +  2 


V„ 


knee 


V^ma x^L  J 


(3-6) 
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where  PDC  is  the  DC  power  dissipated  in  the  device,  Pout  is  AC  output  power  which  is 

R  V2 

given  by  Pout  =— — -  —  —  ,  and  RL  is  a  load  impedance.  For  class  A  operation,  and 

KRon  +Rl) 

reaches  a  maximum  of  50%  for  the  case 

3.  Basic  Dimension  parameters  for  SIT  Design 

We  have  seen  the  SIT  in  Fig.  2. 1  and  Fig.  2.2,  it  is  defined  by  four  geometry  parameters 
(2a,  Lg,  Lsg,  and  Lgd)  and  one  doping  parameter  (ND).  Theoretically,  the  source  gate  distance 

Lsg  should  be  as  small  as  possible  to  minimize  the  source  resistance.  However  Lsg  is  limited  by 

the  requirement  that  a  low  leakage  path  with  sufficient  breakdown  voltage  must  exist  between 
the  gate  and  source.  A  nominal  distance  of  Lsg  =  0.5  \xm  works  well  and  this  value  was  also 

adopted  for  all  simulations.  As  a  side  note  we  mention  that  the  voltage  gain  improves  when  Lsg 
increases  for  fixed  L  ,  because  for  increasing  drain  bias  the  saddle  point  tends  to  move  toward 
and  eventually  reach  the  source,  as  was  shown  in  Fig.  2.7.  An  increased  Lsg  allows  the  saddle 

point  to  move  further  behind  the  gate  and  increase  the  screening  of  the  drain  induced  field  by 
the  gate,  thereby  increasing  the  blocking  gain.  The  trench  width  ( t )  was  fixed  at  one  half  of 
mesa  width  (2a),  and  except  for  a  small  change  in  spreading  resistance  of  the  drain  drift 
region,  it  had  relatively  little  effect  on  the  static  characteristic. 

One  of  the  most  significant  parameters  for  design  is  the  doping  value,  ND.  In  these 
simulations,  we  used  four  different  values,  5 -1015  cm3,  1-1016ct?T3,  5T076  cm3,  and 
2  •  1017cm”3,  spanning  a  range  of  nearly  two  orders  of  magnitude.  Two  lower  doping  values 
correspond  to  devices  used  by  Northrop  Grumman  for  their  pioneering  work  [2]  whereas  the 
two  higher  doping  values  were  used  by  Purdue  University  [3],  The  gate-drain  length  (  Lgd ) 

needs  to  be  tailored  to  each  doping,  and  is  determined  by  the  two  conflicting  requirements 
of  low  drain  resistance  and  a  large  breakdown  voltage.  In  lower  doping  levels,  the  gate- 
drain  length  ( Lgd )  is  limited  by  the  requirements  for  a  low  drain  resistance,  due  to  the  low 

doping  and  resulting  in  high  bulk  resistivity.  On  other  hand,  in  high  doping,  the  bulk 
resistivity  is  low  enough  that  the  gate-drain  distance  can  be  chosen  to  be  comparable  to  the 
punch-through  depletion  depth  for  maximum  drain  voltage.  We  investigate  the  effect  of 
varying  two  design  parameters,  gate  length  ( Lg )  and  mesa  width  (2a)  on  the  static  I-V 

curve  by  keeping  all  other  parameter  constant.  The  computed  static  I-Y  curve  for  SIT  with 
a  doping  value  of  Nd  =  5T075  cm'3  are  shown  Fig.  3.2  for  five  different  gate  lengths  (0.05, 
0.1,  0.2,  0.3,  and  0.5  \xm)  and  five  different  mesa  width  (0.4,  0.5,  0.6,  0.7,  and  0.8  \xm). 
Other  plots  for  different  doping  values  can  be  found  in  appendix  C.  For  the  case  of  a  large 
mesa  width  and  very  short  gate,  as  shown  in  the  lower  left  hand  corner  of  Fig.  3.2,  it  is 
difficult  to  pinch  off  the  drain  current  at  high  gate  voltages,  which  corresponds  to  a  low 
voltage  gain  p.  The  other  extreme  case  of  a  long  gate  length  and  a  narrow  mesa  width,  is 
the  device  characteristic  in  the  upper  right  hand  corner  of  Fig. 3. 2.  This  device  is  pinched 
off  for  all  drain  voltage  with  a  relatively  small  gate  voltage  at  the  expense  of  having  only 
a  very  small  maximum  current  7max . 
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Drain  Current  [ niA/mm ] 


Figure  3.2. 1-V  curves  for  various  gate  length  and  mesa  widths  with  No  =  5*1015 cm\  L,y=0.5  pm and  t=2a.  Curve  start 
at  a  gate  voltage  Vg=0  and  progress  in  -2V  steps. 
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Gate  Length  |/im] 


Fig.  3.5  shows  that  maximum  current,  /max ,  depends  on  the  gate  length,  Lg,  along  with  mesa 
width  2a  as  a  parameter.  In  Fig.  3.5,  the  effect  of  an  increased  Lg  is  a  reduction  in  7max ,  as  a 

result  of  the  increased  channel  resistance.  The  part  of  the  channel  between  the  gates  that  is 
partially  depleted  and  which  forms  the  bottleneck  for  current  flow  simply  becomes  longer.  The 
increase  of  the  channel  resistance  follows  nearly  linearly  with  Lg ,  and  we  can  see  an  inverse 
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proportional  effect  on  the  maximum  current.  On  other  hand,  maximum  current  is  proportional  to 
mesa  width  and  it  increases  with  mesa  width  in  a  nearly  linear  fashion,  because  the  bottleneck 
between  gate  electrodes  becomes  narrower.  The  blocking  gain,  p,  is  a  stronger  function  of  mesa 
width  and  gate  length  than  the  maximum  current,  unless  the  mesa  width  becomes  very  narrow, 
as  shown  Fig.  3.6.  The  physical  reason  is  that  the  current  changes  linearly  with  channel 
dimension  whereas  the  blocking  gain  to  first  order  depends  exponentially  on  the  channel 
dimension  [5],  because  of  high  voltages  and  low  current  in  exponentially  related  to  the  saddle 
point  potential.  As  a  result,  the  blocking  gain  increases  and  the  maximum  current  decreases, 
when  the  ratio  between  gate  length  and  mesa  width  is  increased. 


Gate  length  |/tm| 


Mesa  width  |pm| 


Figure  3.4  Maximum  Current  for  various  gate  lengths  and  mesa  widths. 
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Figure  3.5  Blocking  gain  as  a  function  of  gate  lengths  and  mesa  widths. 

4.  Drain  and  Source  Resistance 

The  parasitic  resistances  in  the  SIT  are  the  resistance  of  the  source  (Rs),  Drain  (RD) 

and  the  bulk  resistance  of  the  substrate  (R™6).  The  substrate  resistance  can  only  be 
meaningfully  defined  for  a  complete  device,  and  by  assuming  current  spreading  angle  of 
45  degrees  can  be  computed  by 


^ sub  P sub 


1 

,(T„t+W)Z 

\w-z\ 

(T,ul+Z)W 

(3.7) 


which  is  similar  to  the  equation  for  the  thermal  resistance  (see  Eqn.  (2.5)),  expect  that  now 


the  resistivity  of  the  substrate,  P sub ,  appears.  The  resistance  of  the  drain  and  source  contact 
for  a  complete  device  can  be  estimated  by 


,  = _ Pd_ _ 

D  (»'+2-r-Xz+2.r„4) 


(3.8) 


R 


s  ~ 


Pd 

2a-W-Nm 


(3.9) 


24 


Typical  values  are  0/m  ^  psub -\5mCl- cm ^  an(j  pDS  cm  ^hen  we 

assume  a  small  area  test  device  with  1  mm  periphery  with  20  fingers,  Z=40  pm  and  W=50 
\im,  the  substrate  contribution  would  be  Rsub  =  2.9Q  •  mm  by  Eqn.  (3.7),  nearly  order  of 

magnitude  smaller  than  the  large  area  limit.  Similarly,  we  find  RD  =  0.5Q  •  mm  and 
Rs  =  0.0024Q  •  cm  for  the  large  area  and  small  area  device  respectively.  From  this  follows 
immediately  that  the  resistance  due  to  the  drain  contact  can  be  neglected;  this  would  still  be  true 
for  even  a  poor  contact  of  /?  =  1  •  1 0  4  Q  •  cm 2 .  The  source  resistance  depends  on  the  particular 
design  because  the  area  of  the  contact  is  determined  by  the  mesa  width. 


5.  Basic  considerations  of  Power  and  frequency  performance 

In  previous  definition  of  maximum  current  7max,  a  large  maximum  current  can  be 
obtained  for  larger  channel  dimensions.  Shown  in  Fig.  3.7  is  the  effect  of  varying  gate 
length  on  frequency  performance.  It  should  be  understood  that  fT  is  a  function  of  zero 
bias.  The  use  of  smaller  dimensions  allows  for  higher  frequency  performance,  but  comes  at 
the  expense  of  a  slightly  lower  maximum  current.  Therefore,  to  maximize  frequency 
performance,  a  short  gate  length  should  be  selected. 


Figure  3.6  Frequency  as  function  of  gate  lengths  and  mesa  widths 
for  fixed  voltage  gain  10  dB 
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CHAPTER  IV 

HIGH  FREQUENCY  CHARATERISTICS 

I.  Small  Signal  Equivalent  Circuit  of  SIT 

The  SIT  is  usually  used  as  a  small  signal  amplifier.  Small  signal  analysis  is  as 
important  as  DC  analysis  for  semiconductor  device,  especially  for  high  frequency  operation 
of  the  SIT.  There  are  three  types  of  techniques  for  small  signal  analysis  of  a  semiconductor 
device  [1],  The  first  technique  is  transient  excitation  followed  by  Fourier  decomposition, 
second  is  incremental  charge  partitioning,  and  third  is  small-signal  AC  steady  state 
frequency  domain  analysis.  The  small-signal  AC  steady  state  frequency  domain  analysis  is  a 
superior  approach  comparing  to  the  other  methods  [1].  This  technique  is  easy  to  be 
incorporated  into  a  device  simulator.  In  this  chapter,  the  small-signal  AC  steady  state 
frequency  domain  analysis  technique  is  applied  on  the  frequency  performance  of  the  SIT. 

For  an  assessment  of  the  frequency  performance  of  the  SIT,  the  small  signals  Y- 
parameters  have  been  computed  by  MEDICI  over  the  anticipated  operation  frequency  range 
from  500  MHz  to  10  GHz.  The  Y-parameters  as  a  function  of  frequency,  as  obtained  from 
MEDICI  simulation,  were  fit  to  the  pi-equivalent  circuit  shown  in  Fig.  4.1.  The  gate-drain 
capacitance  is  occasionally  referred  to  as  the  feedback  capacitance  because  it  couples  the  output 
and  input  of  the  device.  The  model  the  frequency  behavior  of  the  input  impedance  the  gate- 
source  resistance  Rgs  is  connected  in  series  with  Cgs .  The  gate-source  resistance  represents  the 

resistance  of  the  undepleted  part  of  the  mesa  on  the  source  side  of  the  gate. 

The  magnitude  of  the  transconductance,  gm ,  is  equal  to  the  DC  transconductance. 
The  transit  time,  z  ,  accounts  for  the  delay  due  to  the  finite  time  that  a  carrier  needs  to  travel 
from  the  source  to  the  drain.  The  output  conductance  is  denoted  by  Gds .  The  procedure  of 

finding  the  equivalent  circuit  elements  from  the  simulated  data  is  identical  to  determining 
these  elements  from  a  measured  data  set,  with  the  exception  that  measured  data  are  obscured 
by  unwanted  effects  such  as  noise,  imperfect  calibration  and  unaccounted  parasitics  because 
the  exact  details  of  a  physical  structure  are  only  known  within  certain  limits,  whereas  the 
simulated  structure  is  known  precisely.  At  circuit  level,  input  current  7(  is  equal  to  the 

magnitude  of  ideal  output  current  gm  •  Vgs ,  of  the  intrinsic  transistor.  When  the  output  is 
short-circuit  [2], 

i,=Ac„+csiyg,  (4-D 


if  let  Cg  =  Cgs  +  Cgd ,  then  at  the  cutoff  frequency 


71  =  2  *frC,V„  =g„Vil 

(4-2) 

in  Eqn.  (4-2), 

r  §  m 

T  " 

(4-3) 

/r“2dcg,  +  cgJ 

(4-4) 
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The  four  elements  of  the  general  pi-circuit,  7 ,  Y0,  Yr  and  Yf ,  can  be  obtained 

easily  from  the  Y-parameters  [3]  that  one  has  available  either  from  simulations  or  from 
measured  S-parameters  that  were  transformed  to  Y-parameters  (see  App.  D). 


J(.=7ii-7,2  (4-5) 

Y0=Y22-Yu  (4-6) 

Yr=~Yn  (4-7) 

Yf=Y2[-Yn  (4-8) 


Gate  Cgd  Rd  [>rain 


Figure  4.1  Pi-equivalent  circuit  for  small  signal  modeling. 

From  these,  the  reactive  components  of  the  intrinsic  device  may  be  determined.  MEDICI  will 
be  able  to  create  the  Y-parameters  directly  from  the  small  signal  simulation.  Using  circuit 
model  for  the  SIT,  the  value  of  the  components  can  be  determined  by  the  simulations. 
Applying  the  small  signal  simulations  to  the  circuit  models,  from  Eqn.  (4-5)  through  (4-8) 
the  Y-parameters  are  related  by  [3]: 
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Source 

Figure  4.2  General  Hybrid-Pi  equivalent  circuit. 


1)1  =  MCgs+cvl) 

(4-9) 

Y\2  =  ~j0}  ‘  Cgd 

(4-10) 

Y»  =  gm-<rJm-MCgd 

(4-11) 

Y22  =glk  +  jMC*  +  Cgd) 

(4-12) 

The  Y-parameters  used  for  the  following  discussion  are  for  a  bias  point  with  Vd  =10V  and  Vg  — 
0.  Plotting  the  magnitude  of  }/as  a  function  of  frequency,  as  shown  in  Fig.  4.3,  we  can  obtain 
values  for  gm,  and  r .  By  Eqn.  (4-11),  (4-10),  and  (4-8),  the  magnitude  of  |f}|is  almost 

independent  of  frequency,  and  we  obtain  gm  by  forming  the  average  over  the  lower 

frequency  part  of  the  plot.  The  phase  plot  linearly  decreases  over  frequency,  and  r  can 
be  obtained  from  the  slope  of  the  line  by  fitting  the  plot  using  the  equation 
A Yf  =2 7$t  .  Also  shown  in  Fig.  4.3  is  the  magnitude  of  70  that  Gate-Drain  and  Gate- 

Source  Capacitance  from  which  we  obtain  a  value  for  gds.  The  feedback  capacitance 
Cgd  is  obtained  by  averaging  the  data  of  a  plot  of  3(-712  /&>)over  frequency  by  Eqn. 
(4-10),  which  should  be  independent  of  frequency.  Here  3(-712/ry)  means  that  the 
imaginary  part  is  taken  of  -Ynlco  expression  in  brackets.  Similarly,  Cgs  is  found  by 
averaging  3((7n  +  712)/&>)over  frequency  by  Eqn.  (4-9),  (4-10),  as  shown  in  Fig.  4.4. 
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Figure  4.4  Gate-Drain  and  Gate-Source  Capacitance  as  function  of  frequency  for 
bias  point  Vd=  10K,  Vg  =  OK. 

If  one  tried  to  find  a  value  for  Cds  by  simulation  from  MEDICI,  the  resulting  values  seem 
randomly  jump  between  positive  and  negative  numbers  with  an  absolute  value  that  is  at  least 
an  order  of  magnitude  smaller  than  either  of  the  values  found  for  Cgs  or  Cgd .  Therefore,  let 

Cgd=  0.  This  is  physically  meaningful  if  one  considers  the  vertical  arrangement  of  the 
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device  that  essentially  places  the  gate  electrode  as  a  grid  between  the  drain  and  source 
electrode,  electrostatically  shielding  the  two  from  each  other.  This  also  agrees  with  the 
usual  interpretation  of  the  drain  source  capacitance  in  lateral  devices,  such  as  MESFET, 
where  it  is  due  to  fringing  effects  between  the  source  and  drain  electrode  through  the 
substrate  under  the  active  device  [4],  There  is  no  equivalence  to  the  substrate  fringing 
effect  in  a  SIT  due  to  the  vertical  channel,  which  is  similar  to  two  MESFET  channel  in  a 
mirror  image  like  arrangement. 

2.  Effect  of  contact  Resistance 

The  contract  resistivity  of  Rs  and  drain  contact  Rd  can  be  described  by  lumped 
resistances  that  are  in  series  with  the  source  and  drain  terminal,  respectively.  A  finite 
source  resistance  Rs  affects  the  transconductance  as  well  as  the  short  circuit  input 

capacitance.  For  the  special  case  that  Rd=  0  we  find  for  the  transconductance  [5] 

when  Rd=0  (4-13) 

l  +  gmRs 


and  for  the  short  circuit  input  capacitance 


C,„  *  C 


1 


gs 


1  +  SmRs 


■  +  c 


gd 


when  =0 


(4-14) 


The  drain  resistance,  Rd  acts  like  a  built  in  load  resistance  and  therefore  increase  the  input 
capacitance  according  to  the  Miller  effect.  When  we  assume  in  any  special  case  that 
R  =  0 

s  ,  the  Miller  effect  can  be  theoretically  expressed  by  the  following  approximation  [5] 


r 


c  ~c  +c , 

in  gs  ^  gd 


Sn 


\ 


V  §ds+l/RdJ 


assumed  R  =  0 


(4-15) 


In  the  general  case,  when  both  Rs  and  Rd  are  unequal  to  zero,  the  expression  for  the  input 


capacitance  becomes  more  complicated.  However,  we  conclude  that  an  increase  in  Rd  and  Rs 
leads  to  an  overall  degradation  in  frequency  performance,  namely  the  cut-off  frequency,  because 
the  contribution  of  Cgd  is  significant  for  the  SIT  and  tends  not  to  decrease  when  a  source 

resistance  is  introduced,  Eqn.  (4-14).  Note  that  for  a  device  with  a  very  small  feedback 
capacitance,  the  cut-off  frequency  does  not  change  notably  with  increasing  source  resistance. 
Such  as,  MESFET  and  HEMT,  where  the  overall  contribution  of  Cgd  to  the  total  input 

capacitance  is  typically  less  than  10%  [6].  The  source  resistance  is  also  increases  the  total 
resistance  of  the  device  between  source  and  drain,  and  therefore  the  maximum  current  reduced. 
The  effect  of  the  specific  contact  resistivity  pc  is  illustrated  from  Fig.  4.5  to  Fig.  4.15. 
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Drain  Resistivity  =  0.5*10'5  [ft  cm2] 


Figure  4.5  Cut-off  Frequency  as  a  function  of  gate  length  and  mesa 
drain  contact  resistivity  0.5-10'5  |Q  cm2]. 
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Drain  Resistivity  =  1*1 0~6  [£}  cm2] 
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Figure  4.6  Cut-off  Frequency  as  a  function  of  gate  length  and  mesa  width 
with  drain  contact  resistivity  1*1(T5  [fl  cm2]. 


Drain  Resistivity  =  6*1  O'*  [O  cm2] 
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Figure  4.7  Cut-off  Frequency  as  a  function  of  gate  length  and  mesa  width 
with  drain  contact  resistivity  5*105  [O  cm2|. 
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Drain  Resistivity  =  10*1  O'8  [li  om2) 
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Figure  4.8  Cut-off  Frequency  as  a  function  of  gate  length  and  mesa  width  with 
drain  contact  resistivity  1(M0'S  |ft  cm2]. 

Drain  Resistivity  =  60*10*6[n  cm2] 
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Figure  4.9  Cut-off  Frequency  as  a  function  of  gate  length  and  mesa  width  with 
drain  contact  resistivity  50*10"s  |Q  cm2). 
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Figure 


Sourot  Resistivity  =  1e-005[u  cm2] 
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4.10  Gate-Drain  Capacitance  as  a  function  of  gate  length  and  mesa  width  at 
source  contact  resistivity  1-10"5  [ft  cm*],  from  top  drain  contact 
resistivity  50*  KT5, 10*10“5,  5*10'5, 1-10"5, 0.5-105  [ft  cm2],  respectively. 


Source  Resistivity  =  1e-006[Q  cm*] 
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Figure  4.11  Gate-source  Capacitance  as  a  function  of  gate  length  and  mesa  width 
at  source  contact  resistivity  1-10"5  [ft  cm2),  from  top  drain  contact 
resistivity  5(H0‘S,  10-105,  5*1  O'5,  MO  ,  0.5*1 0  s  jfl  cm2|,  respectively. 
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Source  Resistivity  =  1o-006[u  cm2] 


Figure  4.12  Transconductance  as  a  function  of  gate  length  and  mesa  width  at 
source  contact  resistivity  MO'5  |il  cmz|,  from  top  drain  contact 
resistivity  0.5- 10  5,  1  •  1  O'5,  5*10's,  1 0- 1 O'5,  50*10's  [ft  cm2),  respectively. 
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Figure  4*13  Cut-off  Frequency  as  a  function  of  gate  length  and  mesa  width  at 
source  contact  resistivity  1*1 0“5  |0  cm2],  from  top  drain  contact 
resistivity  0,5*10“5, 1-10"5,  5-10"5, ltMO5, 50- 10*5  [Q  cm2],  respectively. 
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Drain  Resistivity  =  1*1  O’6  [£i  cm2] 
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Figure  4.14  Transconductance  as  a  function  of  gate  length  and  mesa  width  at  drain 
contact  resistivity  1*10'5  [ft  cm2],  source  contact  resistivity  0.5*10'5, 1*1(T, 
5*10‘5, 10-10'5,  50*1  O’5  |ft  cm2],  from  top  to  bottom,  respectively. 
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Figure  4.15  Gate-drain  capacitance  as  a  function  of  gate  length  and  mesa  width  at 
drain  contact  resistivity  1*10"5  [ft*cm2|,  source  contact  resistivity  50*10’5, 
10*10‘5,  5-1  O’5,  MO’5,  0.5*1 0‘5  [ft  cm2|,  from  top  to  bottom,  respectively. 
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Drain  Resistivity  =  1*10'5  [fi  cm2] 
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Figure  4.16  Gate-source  capacitance  as  a  function  of  gate  length  and  mesa  width  at 
drain  contact  resistivity  1*10'5  [Q*cm2|,  source  contact  resistivity  0.5*10  , 
MO'5, 5*10'5, 10-10'5,  50- 10  s  (Q  cm2J,  from  top  to  bottom,  respectively. 


3.  Effect  of  Mobility 


The  carrier  mobilities  and  account  for  scattering  mechanisms  in  electrical  transport.  MEDICI 
provides  several  mobility  model  choices. 

As  an  alternative  to  the  concentration-dependent  mobility,  this  mobility  is  concentration- 
and  temperature-dependent  empirical  mobility  function.  These  are  given  by  the  expressions  [7] 


f  T 

f^mtx  I  200  J  7^ mm 

M"P  =  +  U(T 

300;  N^p 


(4-16) 


In  most  case,  this  mobility  model  use  in  semiconductor  device  simulation. 

3.1  Perpendicular  Electric  Mobility  Model 

This  model  is  applied  at  every  position  in  the  device  and  not  just  at  interfaces  or  in 
the  inversion  layer  and  is  described  by  modified  version  of  Analytic  Mobility  Eqn.  (4-16). 
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(4-17) 


where  1 " ,  and  L’p  are  the  components  of  electric  that  are  perpendicular  to  the  side  of  an 
element  (the  default)  or  the  components  of  electric  perpendicular  to  the  current  direction. 
G  G 

The  factors  surfacn  and  surfacp  are  only  applied  at  interfaces  between  semiconductor  and 
insulator. 


3.2  Surface  Mobility  Model 

Along  insulator-semiconductor  interface,  the  carrier  mobilities  can  be  substantially 
lower  than  in  the  bulk  of  the  semiconductor  due  to  surface  scattering.  This  model  is  only 
applied  at  the  interface  because  it  assumes  that  the  carrier  inversion  layer  width  is  smaller 
than  the  grid  spacing  used  at  the  interface.  With  this  assumption,  all  the  inversion  charge 
effectively  occurs  at  the  interface,  and  it  is  appropriate  to  use  this  model.  This  model 
calculates  effective  mobilities  at  the  interfaces  using  the  expressions 


Ms,  n,p  (Surf ace)  ^ surf an, p 


Jeffl,n,p 


(4-18) 


Also,  and  enhanced  surface  mobility  model  has  been  included  in  MEDICI  that  takes  into 
account  phonon  scattering,  surface  roughness  scattering,  and  charges  impurity  scattering. 
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The  Hewlett-Packard  Mobility  model  takes  account  dependence  on  electric  s  both 
parallel  and  perpendicular  to  the  direction  of  current  flow.  The  expressions  for  mobility 
used  by  this  model  are 


M„, 


Ml n.f 


1  + 


f  p.  \2 

"  y"p 

\ _ <■"  _ y 

ML,n,pEujt,p  i  /-HP 

yHP  ' 


(4-19) 


f  Ml.n.p^O.n.i 
yHP 
-m.p  y 


The  expressions  for  juXnp  are  given  by 


Mi  ,,.P  = 


 Mo 


HP 


E,„ 

1  +  ^hP 


if  Nmai(x,y)  <  N% 


E 


HP 

cnyp 


(4-20) 


the  default  value  for  N"£  is  5-1017.  if  above  condition  is  not  satisfied,  then  m±.„,p  =  M,hp , 
where  fin  p  are  low  mobility  values  described  in  4.3.1  Analytic  Mobility. 
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Figure  4.17  Cut-off  Frequency  depend  on  mobility  model 
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Figure  4.18  Cut-off  Frequency  in  several  mobility  models 
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4.  Effects  of  Saturation  Drift  velocity 

In  semiconductor  device  design,  one  of  the  most  important  parameters  is  a  saturated 
drift  velocity  that  the  relationship  between  the  drift  velocity  of  electrons  and  applied 
electric  field.  The  first  measurements  of  the  high-field  drift  velocity  of  electrons  in  SiC 
were  performed  by  von  Muench  and  Pettenpaul  in  1977  on  Lely  crystals  of  6H-SiC.  They 
reported  an  electron  saturation  velocity  of  2x10 1  cm/s  parallel  to  the  basal  plane  at  room 
temperature.  But  recently,  the  drift  saturation  velocity  reported  to  1.9x10 1  cm! s  in  6H-SiC 
and  2.2x10 1  cml s  in  4H-SiC,  with  data  at  both  room  temperatures  [8].  Therefore,  in  this 
work,  we  simulated  device  behavior  with  different  saturated  drift  velocity,  respectively. 
Shown  in  Fig.  4.19. 
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Figure  4.19  Cut-off  Frequency  depend  on  Saturation  Drift  velocity. 
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Figure  4.20  Transconductors  depend  on  Saturation  Drift  velocity,  from  top 
Vjut  is  2.5  x  107,  2.2  x  10%  2.0  x  107  and  1.5  x  10 1\cm/s\. 
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Figure  4.21  Gate-drain  capacitances  in  other  saturation  drift  velocity,  from  top  vsat 
is  2.5  x  107, 2.2  x  107, 2.0  x  107,  and  1.5  x  10 7(cm/s]. 
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5.  Considerations  for  higher  frequency  SIT  Design 

In  order  to  design  SITs  with  higher  frequency,  scaling  of  the  device  doping  and  mesa  with 
has  the  biggest  effect.  The  transconductance  increase  faster  with  doping,  roughly  linearly,  than 
the  input  capacitance.  The  cut-off  frequencies  as  a  function  of  doping  and  for  gate  length, 
assuming  ideal  contacts,  shown  in  Fig  4.22,  cut-off  frequencies  in  excess  of  20  GHz  can  be 
obtained  for  doping  2-10 1  cm^ .  However,  reduction  of  trench  width  is  also  important  because 
the  value  of  Cgd  _  which  contributes  to  the  input  capacitance  Cm  =  Cgd  +  Cgs  of  device, 

essentially  scales  with  the  trench  width  t.  The  effect  of  Cgd  gains  its  importance  when  the 
doping  is  increased  because  Cgd  can  be  approximated  by  the  depletion  capacitance  under  the 

gate  electrode.  Also,  the  reduction  in  depletion  width  with  increasing  doping  for  a  given  bias 
point  results  in  an  increase  value  for  Cgd .  For  a  scaled  trench,  the  effect  of  Cgd  eventually 

becomes  so  pronounced  that  the  rise  in  cut-off  frequency  for  increasing  doping  levels  off. 
Therefore,  it  is  important  to  scale  the  trench  together  with  the  mesa  width.  However,  other 
considerations  besides  maximizing  the  frequency  response  affect  the  choice  of  the  trench  width, 
such  as  device  heating  due  to  excessive  power  density. 


Figure  4.22  Cut-off  Frequency  as  a  function  Doping  No  and  Gate  length. 
Specific  contact  resistance  for  source  is  zero. 
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CHAPTER  V 

TEMPERATURE  SIMULATIONS 

1.  Theoretical  Considerations 

Silicon  Carbide  (SiC)  is  a  semiconductor  material  combining  several  unique 
features  such  as  high  breakdown  field,  high  electron  saturation  velocity,  and  high 
thermal  conductivity.  These  properties  make  it  a  promising  candidate  for  RF  power 
transistors.  RF  Power  transistors  are  commonly  operated  at  high  voltage  and  high 
current  levels  leading  to  considerable  self-heating.  Thus,  the  operating  temperature 
can  be  much  higher  than  room  temperature.  To  simulate  power  transistors  under 
these  conditions  properly,  the  temperature  dependence  of  the  electron  mobility  in 
the  both  the  low-field  and  high-field  region  must  be  also  taken  into  account  in  the 
mobility  models.  Recently,  the  electron  mobility  in  4H,  6H,  and,  3C  SiC  has  been 
investigated  both  experimentally  and  theoretically  by  Monte  Carlo(MC)  simulations 
by  several  authors,  and  models  for  the  electron  mobility  have  been  proposed. 
Unfortunately,  none  of  these  models  includes  detailed  temperature  dependence  for 
both  the  low-field  and  high-field  mobilities  [1].  In  this  works,  we  will  use  the 
electron  mobility  model  for  SiC  RF  device  simulation.  The  models  are  based  on 
reported  mobility  data  and  include  the  dependence  of  the  mobility  on  doping 
concentration,  temperature,  and  electric  field.  The  following  4H-SiC  temperature 
dependant  parameters  have  been  determined  by  [1] 
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2.  Gate-Drain  Breakdown  Voltage 

Silicon  Carbide  (SiC)  has  received  increasing  attention  for  power  switching, 
microwave  and  high  temperature  applications  due  to  its  high  breakdown  electric  field, 
thermal  conductivity  and  electron  saturation  drift  velocity.  One  of  the  most  important 
parameters  of  SiC  power  device  is  its  breakdown  voltage.  Also,  one  of  the  advantages  of 
the  static  induction  transistor  is  its  high  output  impedance  which  results  from  the  high 
breakdown  voltage  between  the  gate  and  drain  terminals.  This  is  expected  because  of  the 
vertical  arrangement  of  the  device  channels.  The  gate  electrode,  in  effect,  looks  like  a 
finely  meshed  grid  such  that  the  overall  field  distribution  approaches  that  of  a  parallel 
plate  capacitor.  The  local  field  at  the  very  comer  of  the  gate  is  expected  to  approach 
arbitrarily  high  values,  depending  on  how  sharp  the  comer  will  be  in  practice,  as  shown 
in  Fig.  5.1  through  Fig.  5.3.  However,  this  field  spike  does  not  mean  that  the  breakdown 
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due  to  avalanche  will  occur  at  much  lower  voltages  than  predicted  by  calculation.  The 
avalanche  breakdown  process  occurs  when  electrons  and/or  holes,  moving  across  the 
space  charge  region,  acquire  sufficient  energy  from  the  electric  field  to  create  electron- 
hole  pairs  by  colliding  with  atomic  electrons  within  the  depletion  region.  Thus,  for 
avalanche  to  occur,  a  high  enough  electrical  field  must  be  supplied  to  the  charge  carriers 
in  order  to  produce  electron-hole  pairs,  and  a  certain  extent  of  the  avalanche  region  is 
needed  to  obtain  a  round  trip  gain  of  unity  in  order  to  sustain  the  avalanche.  In  this  work, 
because  of  the  importance  of  the  gate-drain  breakdown  voltage,  the  maximum  obtainable 
output  voltage  for  a  SIT  was  investigated  by  means  of  simulations  with  several 
temperatures  and  doping,  respectively. 


Figure  5.1  Avalanche  Breakdown  voltage  and  Electric  Field  at  Impact  Ionization 
integral  reaches  unity  for  300K  and  400K. 
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Figure  5.2  Avalanche  Breakdown  voltage  and  Electric  Field  at  Impact  Ionization 
integral  reaches  unity  for  500K  and  600K. 
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Figure  5.3  Avalanche  Breakdown  voltage  and  Electric  Field  at  Impact  Ionization 
integral  reaches  unity  for  700K  and  800K. 


In  Fig.  5.1  to  Fig.  5.3,  we  can  see  that  avalanche  breakdown  voltages  decrease, 
maximum  electric  field  increase,  and  depletion  region  increase.  The  breakdown  voltage 
due  to  avalanche  breakdown  predicted  by  the  ionization  integral  is  usually  lower  for  the 
SIT  because  it  does  not  take  into  account  that  a  carrier  can  escape  from  the  path  along  the 
highest  field  either  by  diffusion  or  scattering  due  to  the  curvature  of  the  electric  field  line, 
thereby  decrease  the  multiplication  factor  of  carriers  in  the  avalanche  region.  In  Fig.  5.4 
through  Fig.  5.7,  four  different  drain  regions  are  considered  with  doping  concentrations  of 
5-1015cm~3,  l-1016cm~3,  5 -1016  cm3,  and  2 -1017  cm3  that  have  a  thickness  of  gate-drain, 
Lgd=4[im,  3 pm,  2pm,  and  1.3pm,  respectively.  When  the  gate-drain  voltage  is  not  limited 

by  the  gate  edge,  the  device  will  break  down  intrinsically  at  the  corner  of  the  gates  facing 
the  mesa  center.  The  simulated  breakdown  due  to  avalanche  of  the  device  668V,  473V, 
382V,  and  174V  for  the  two  considered  drain  regions,  and  the  avalanche  breakdown 
voltage  is  roughly  25%  lower  than  one  would  obtain  from  a  parallel  plate  calculation. 
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Figure  5.4  Avalanche  Breakdown  voltage  with  temperature  increasing  for  the  doping  5-1016  \cni3\ 
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Figure  5.5  Avalanche  Breakdown  voltage  with  temperature  increasing  for  the  doping  2*1017  [cm3\ 
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Figure  5.6  Avalanche  Breakdown  voltage  with  temperature  increasing  for  the 
doping  5-1016  [cm3 j,  from  top  300K,  400K,  500K  600K  700K,  800. 
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Figure  5.7  Avalanche  Breakdown  voltage  with  temperature  increasing  for  the 
doping  2*1017  \cni3\,  from  top  to  bottom  300K,  400K,  500K  600K  700K, 
800K. 
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3.  Temperature  Dependent  DC  performances 

DC  performances  at  the  room  temperature  have  been  described  in  chapter  3.  Fig  5.6 
shows  the  I-V  curves  from  300K  to  800K.  As  the  operating  temperature  increases,  Ids 

decrease  monotonically  for  Vg  =  OF.  At  700K,  the  saturation  current  it  almost  50%  of  the 

300K.  This  is  quite  good  comparison  with  Si  device,  which  is  completely  useless  at  such 
temperature.  The  temperature  dependent  characteristics  can  be  explained  as  follows.  As 
temperatures  increase,  more  donors  are  ionized;  the  depletion  width  decreases,  and  the 
effective  channel  thickness  increase.  Thus,  the  total  number  of  mobile  carriers  increases. 
On  the  other  hand,  the  carrier  mobilitizes  decreases.  The  maximum  current  is  shown  in 
Fig.  5.7.  As  temperature  increases,  the  maximum  current  decreases. 
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Figure  5,8  I-V  Curves  at  (a)  T=300tf,  (b)  T=400A,  (c)  T=500A\  (d)  T=600tf,  (e) 
T=700A,  (f)  T=800if,  respectively.  From  the  top  to  the  bottom,  the  five 
curves  in  each  figure  correspond  to  Vg=0K  -2F,  -4F,  -6F,  -8F,  -lOP,  - 
12K-14K-16K-18K 
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Figure  5.9  I-V  Curves  at  T=300K,  400K,  500K,  600K,  700K,  800K  for 
Vg=8V,  10  V,  from  the  top  to  the  bottom,  respectively. 
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igure  5.10  Maximum  Current  at  T=300iT,  400JC,  500X,  600/f,  700A",  800A  for 
2'10l7cm3,  from  the  top  to  the  bottom,  respectively. 
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4.  Temperature  Dependent  Frequency  performances 

Since  the  system  input  and  output  impedance  are  variables,  we  are  interested  in 
the  maximum  frequency  that  the  device  can  operate  The  definition  of  this  maximum 
frequency  becomes  ambiguous.  These  are  two  standard  measures  of  this  maximum 
frequency;  fT  and  /max .  The  fT ,  or  the  frequency  bandwidth  is  the  frequency  that  the 
transistor  has  unity  current  gain  while  driving  a  short  circuit  load.  This  is  equal  to  the 
h2X  =  1  frequency.  The  /max  is  the  maximum  frequency  of  oscillation.  In  order  to  have 

oscillation,  the  gain  of  the  device  must  be  equal  to,  or  above  unity.  Therefore  /max  is  the 
highest  frequency  unity  gain  can  be  obtained  for  any  load  condition.  This  is  obtained 
under  the  maximum  available  gain  condition;  ( Gmax )  which  occur  when  the  device  is 
conjugate  matched. 

Y  or  5-parameters  can  be  extracted  from  MEDICI  simulation  and  h2l ,  Gmax  can 
be  obtained  by  calculation.  The  h2X  can  be  calculated  by 


h2l  — 


■25 


21 


1  21 
4 


(5-6) 


This  is  using  standard  parameter  transformations.  From  this  fT  is  obtained.  The 
maximum  available  power  gain  can  also  be  calculated  from  Y  or  5-parameters. 


G 


max 


where  K  is  the  stability  factor,  which  is  given  by 


K  = 


^11^21  ^12^21 

2 

|S„ 

2-|s!; 

2 

2|S,2| 

52[ 

(5-7) 


(5-8) 


The  stability  factor  will  predict  the  stability  of  the  device  against  oscillation.  For 
a  K  greater  than  unity  and  positive,  the  device  will  be  unconditionally  stable  for  any 
impedance  presented  to  it.  Oscillations  are  possible  whenever  an  |5n|or  |522|  greater 

than  unity.  The  device  is  unconditionally  stable  if  the  input  and  output  impedances  are 
greater  than  zero  for  any  source  or  load  impedance  presented.  The  calculation  of  Gmax  is 
only  valid  when  the  device  is  unconditionally  stable.  From  Eqn.  (5-6)  the  gain  bandwidth 
fT  can  be  calculated  directly  from  the  simulated  Y-parameters.  Fig.  5.8  shows  h2l 
verses  frequency.  The  extrapolated  intercept  of  h2l  =  unity  (0  dB)  can  be  found  fT .  For 
this  SIT,  0.5  \xm  mesa  width  and  gate  length  0.3  \xm  with  Vg  =  -4V ,  Vd  =  50V ,  fT  is 
9.6  GHz  at  600  [K], 
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Magnitude  h2i  [dB] 


Doping  5.0e+16  cm'3,  2a  =  0.5 pm,  Gate  =  0.3  ftm 

26 


Frequency  [GHz] 


Figure  5.11  h2i  parameter  from  simulated  Y-parameter  showing  the  bandwidth- 

gain  product  <fT)  of  9.6  GHz  at  600  [K|.  Bias  point  is  Vj  =  50  V,  Vg=  -4  V. 
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Additionally  from  the  Y-parameters  obtained  by  MEDICI  Simulations,  the  data  can  be 
converted  into  S-parameters,  and  using  information  from  Eqn.  (4-9)  through  (4-12),  the 
circuit  model  capacitance  can  be  determined.  Also,  we  can  obtain  cut-off  frequency  from 
the  capacitance.  Fig.  5.10  thorough  Fig.  5.20  show  cut-off  frequency,  transconductance, 
gate-source  capacitance  and  gate-drain  capacitance  for  a  bias  point  of  Vd  =10,  50,  75, 
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Cut-off  Frequency  [GHz] 


[ V\  and  Vg  =0,  -4  [V\  with  doping  value  5.0e  +  16,  2.0e  +  17  [cm  3],  As  temperature 

increases,  frequency  decreases.  Mainly,  transconductance  affect  the  frequency.  The  other 
side,  gate-drain  capacitance  is  an  inversely  proportional  temperature  increase. 


Figure  5.13  Cut-off  frequency  plotted  as  a  function  of  gate  length  and  mesa 

width  for  a  bias  point  of  Vd  =  50  V,  Vg  =-4  V,  Doping  value  5.0e  +  16cm“3 

from  top  to  bottom  curves,  300,  400,  500,  600,  700,  800  [K], 
respectively. 
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Figure  5.14  Transconductance  plotted  as  a  function  of  gate  length  and  mesa 
width  for  a  bias  point  of  V*  =  50  V,  Vg=  -4  V,  Doping  value  5.0e+16 
cm'3,  from  top  to  bottom  curves,  300,  400,  500,  600,  700,  800  [K|, 


Mesa  Width  fom] 


Figure  5.15  Gate-Source  capacitance  plotted  as  a  function  of  gate  length  and 
mesa  width  for  a  bias  point  of  Vj  =  50  V,  Vg  =  -4  V,  Doping  value 
5.0e+I6  cm'3,  from  top  to  bottom  curves,  300,  400,  500,  600,  700,  800 
[K],  respectively. 
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Figure  5.16  Gate-Drain  capacitance  plotted  as  a  function  of  gate  length  and  mesa 
width  for  a  bias  point  of  Vj  —  50  V,  Vg=  -4  F,  Doping  value  S.0e+16 
cm3,  from  top  to  bottom  curves,  800,  700,  600,  500,  400,  300  [KJ, 
respectively. 
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Figure  5.17  Cut-off  frequency  plotted  as  a  function  of  gate  length  and  mesa  width 
for  a  bias  point  of  Vd~ 25  V,  Vg=  0  F,  Doping  value 5.0e+l6  cm3,  from 
top  to  bottom  curves,  300, 400, 500, 600, 700, 800  JK),  respectively. 
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Figure  5.18  Transconductance  plotted  as  a  function  of  gate  length  and  mesa  width 
for  a  bias  point  of  Vd  =  25  V,  Vg=  0  V,  Doping  value  5.0e+16 cm'3,  from 
top  to  bottom  curves,  300,  400,  500,  600,  700,  800  [K],  respectively. 


Gate  length  fom]  0.4  Mesa  Wjdth 

Figure  5.19  Gate-Source  capacitance  plotted  as  a  function  of  gate  length  and 
mesa  width  for  a  bias  point  of  Vd  =  25  V,  Vg  —  0  K,  Doping  value 
5.0e+16  cm\  from  top  to  bottom  curves,  300,  400,  500,  600,  700,  800 
[KJ,  respectively. 
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Figure  5.20Gate-Drain  capacitance  plotted  as  a  function  of  gate  length  and 
mesa  width  for  a  bias  point  of  Vj  =  25  V,  Vg  =  0  V9  Doping  value  5.0e+16 
cm \  from  top  to  bottom  curves,  800,  700,  600,  500,  400,  300  [K], 
resoectivelv. 
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Figure  5.21Cut-off  frequency  plotted  as  a  function  of  gate  length  and  mesa  width 
for  a  bias  point  of  Vj  =  75  V,  Vg=  -10  V,  Doping  value  S.0e+16  cm'3, 
from  top  to  bottom  curves,  300,  400,  500,  600,  700,  800  [K), 
respectively. 


65 


Gate  length  fom]  0.4  Mesa  wjdth 

Figure  5.22  Transconductance  plotted  as  a  function  of  gate  length  and 
mesa  width  for  a  bias  point  of  V4  =  75  V,  Vg  =  -10  V,  Doping 
value  S.0e+16  cm3.,  from  top  to  bottom  curves,  300,  400,  500, 
600, 700, 800  [K],  respectively. 


Figure  5.23  Gate-Source  capacitance  plotted  as  a  function  of  gate  length  and 
mesa  width  for  a  bias  point  of  Vj  =  75  V,  Vg  =  -10  V,  Doping  value 
5.0e+16  cm'3,  from  top  to  bottom  curves,  300,  400,  500,  600,  700,  800 
[K|,  respectively. 
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Figure  5.24  Gate-Drain  capacitance  plotted  as  a  function  of  gate  length  and  mesa 
width  for  a  bias  point  of  Vj  =  75  V,  Vg  =  -10  F,  Doping  value  5.0e+16  cm3, 
from  top  to  bottom  curves,  800,  700, 600,  500,  400,  300  [K],  respectively. 

1 1 4 


0.046 

•  0.04 


0.035 


0.03 

j^j°.026 


Mesa  Width  lan ] 


Gate  length  [jan] 


Figure  5.25  Cut-off  frequency  plotted  as  a  function  of  gate  length  and  mesa 
width  for  a  bias  point  of  Vj=  10  V,  Vg=  0  V,  Doping  value  1.0e+16  cm 
3 ,  from  top  to  bottom  curves,  300,  400,  500,  600,  700,  800  [K], 
respectively. 
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Figure  5.26  Cut-off  frequency  plotted  as  a  function  of  gate  length  and 
mesa  width  for  a  bias  point  of  Vj  =  10  V,  Vg=  0  V. Doping  value 
5.0e+15  cm'3,  from  top  to  bottom  curves,  300,  400,  500,  600,  700, 
800  IKI.  resDectivelv. 
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Figure  5.27  Cut-off  frequency  plotted  as  a  function  of  gate  length  and  mesa  width 
for  a  bias  point  of  Vj  =  10  V,  Vg=  0  V,  Doping  value  2.0e+17 cm'3,  from 
top  to  bottom  curves,  300, 400, 500,  600, 700, 800  [K],  respectively. 
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Figure  5.28  Transconductance  plotted  as  a  function  of  gate  length  and  mesa  width 
for  a  bias  point  of  Vj  =  10  V,  Vg=  0  V,  Doping  value  2.0e+17  cm'3,  from 
top  to  bottom  curves,  300,  400, 500, 600,  700, 800  [K|,  respectively. 
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Figure  5.29  Cut-off  frequency  plotted  as  a  function  of  gate  length  and  mesa  width 
for  a  bias  point  of  Vj  =  10  V,  Vg=  0  V,  Doping  value  2.0e+17  cm'3,  with 
Source  contact  5e-5  \ilcm2\,  from  top  to  bottom  curves,  300,  400,  500, 
600,  700,  800  [K],  respectively. 
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Figure  5.30  Transconductance  plotted  as  a  function  of  gate  length 
and  mesa  width  for  a  bias  point  of  Vj  =  10  V,  Vg=  0  V , 
Doping  value  2.0e+17  cm'3,  with  source  contact,  5e-5 
[flcm2\,  from  top  to  bottom  curves,  300,  400,  500,  600, 
700, 800  [K],  respectively. 


In  above  results,  as  temperature  is  increased,  fr  and  g„,  are  decreased.  At  800  K,  gm  is 
about  20%  of  its  value  at  300  K.  The  monotonic  decrease  of  transconductance  is  due  to 
the  decrease  of  the  saturation  current  which  is  mainly  a  result  of  mobility  degradation.  At 
800  K  with  doping  2.  0e+ 1 7  cm3,  fr  is  about  1 3  GHz  which  is  about  20%  of  its  value  at 
300  K.  Thus,/r  and  gm  decrease  to  almost  same  rate. 
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CHAPTER  VI 

CONCLUSION  AND  FUTURE  WORK 

1.  Conclusion 

The  superior  physical  properties  of  SiC,  such  as  wide  band  gap,  high  saturation  electron 
velocity,  low  impact  ionization  rates,  and  high  thermal  conductivity,  provide  very 
promising  candidates  for  high  temperature,  high  power,  high  frequency,  and  radiation  hard 
applications.  The  high  blocking  voltage,  the  low  harmonics  and  the  low  internal 
modulation  distortion  (IMD)  are  the  excellent  features  of  SIT  compared  to  those  of  BJT 
and  FET.  The  current-voltage  (I-F)  transfer  function  (A to  Ic)  of  BJT  is  the  exponential, 

and  the  square  law  ( Vg  to  Id )  of  the  conventional  FET.  In  contrast  to  those  of  the  BJT  and 
FET,  the  low  distortion  amplifier  was  reported  by  virtue  of  its  linear  V  -  Id  characteristics 

of  SIT.  In  this  work,  the  static  induction  transistor  (SIT)  in  SiC  was  investigated  as  a 
possible  microwave  device  for  high  breakdown  voltage  and  large  bandwidth  of  operation. 
A  large  number  of  device  structures  were  considered  by  means  of  simulations.  The  key 
design  parameters  are  the  mesa  width  2 a  and  Lg  while  maintaining  adequate  blocking  gain. 

It  is  also  shown  that  short  gate  lengths  are  needed  for  a  high  cut-off  frequency  and  that  a 
trade-off  exists  between  frequency  performance  and  current  density. 

The  main  approach  to  get  higher  cut-off  frequency  is  to  increase  the  doping  in  channel, 
and  adjust  the  mesa  width  and  gate  length  accordingly.  Also,  trench  width  is  made  very 
narrow.  For  the  highest  doping  values,  the  smallest  trench  width  is  considered;  a  cut-off 
frequency  of  up  to  42  GHz  is  predicted. 

Actually,  current  high  cut-off  frequency  devices  do  not  reach  theoretical  best 
performance  because  of  the  source  contact  resistivity  and  the  shrink  source  contact  area 
with  reduced  mesa  width.  Contacts  better  than  1.5  5  Clem2  are  needed  to  reduce  the  effect  of 
the  source  resistance.  The  large  gate-drain  capacitance  of  the  SIT  does  not  only  reduce  the 
cut-off  frequency  but  more  importantly,  the  input  capacitance  under  realistic  operating 
conditions  is  further  increased  due  to  the  Miller  effect. 

The  equivalent  circuit  was  extracted  from  simulated  Y-parameter  as  a  function  of  bias, 
and  it  can  be  expanded  to  a  large  signal  equivalent  circuit.  This  circuit  model  was  imported 
into  a  circuit  simulation  program  to  simulate  load-pull  results,  and  could  also  be  used  for 
circuit  design. 

The  drain-gate  breakdown  voltage  is  limited  by  breakdown  at  the  edge  of  the  gate 
metal,  unless  an  edge  termination  is  employed.  The  edge  termination  of  the  gate  seems  to 
be  of  limited  use  because  the  drain-gate  breakdown  merely  shifts  to  the  inner  trench 
corners,  and  the  ultimate  breakdown  voltage  is  only  slightly  higher  than  in  the 
unterminated  device.  However,  it  is  still  useful  to  terminate  the  gate  edge,  especially 
because  the  added  processing  is  minimal  and  easy  to  perform.  The  breakdown  with 
terminated  gate  is  less  abrupt,  because  a  large  leakage  current  signals  that  the  maximum 
operating  voltage  is  approached  and  the  failure  may  be  thermal  an  therefore  time 
dependant.  This  made  devices  with  termination  more  rebust  during  operation  and  testing. 

The  presented  results  demonstrate  that  simulations  are  accurate  enough  to  be  used  as  a 
design  tool,  and  to  establish  a  reference  line  against  which  measured  device  performance 
can  be  compared.  It  needs  to  be  understood  that  the  results  for  absolute  numbers  need  to  be 
interpreted  with  care,  because  the  material  parameters  of  SiC  are  not  known  well  enough. 
They  may  even  vary  from  lot  to  lot,  to  vendor  to  vendor,  particularly  when  being  compared 
with  the  control  available  in  more  mature  material  systems.  A  better  knowledge  base  will 
emerge,  once  the  commercialization  of  SiC  based  devices  progresses  further. 
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The  two  key  advantages  of  the  SIT  are  the  high  breakdown  voltage  and  the  high 
power  density  per  area.  Both  of  these  advantages  allow  the  SIT  with  large  periphery  for 
large  output  power.  However,  with  increasing  doping  the  advantages  of  a  high  breakdown 
voltage  diminishes,  the  large  breakdown  voltage  only  helps  on  the  output  of  the  device. 
The  bandwidth  is  limited  by  the  large  input  capacitance  due  to  Cgd  and  Cgs .  Also,  the 

inability  to  scale  the  device  to  large  periphery  at  high  frequencies  diminishes  the 
advantages  of  a  large  power  per  unit  area.  Therefore,  a  trade-off  exists  between  high 
breakdown  voltage  and  high  frequency. 

The  SIT  is  the  suitable  device  for  the  recent  digital  RF  systems  such  as  nl  4  shifted 
quadric  phase  shift  keying  (QPSK)  and  quadric  amplitude  modulation  (QAM)  and 
broadband  amplifier  such  as  for  the  CATV. 

2.  Future  Work 

The  neural  networks  have  recently  been  introduced  to  the  microwave  area  as  a  fast 
and  flexible  vehicle  to  microwave  modeling,  simulation  and  optimization.  Also,  in  recent 
years,  a  new  CAD  approach  based  on  neural  network  models  has  been  introduced  for 
microwave  device  design.  It  has  been  applied  to  the  efficient  modeling  of  devices,  circuit 
and  device  optimization. 

Even  if  today’s  simulations  are  powerful  tools,  they  can  only  be  as  accurate  as  the 
mathematical  models  that  are  used.  As  a  consequence,  constructing  good  model  for  device 
used  in  circuits  for  circuits  used  in  systems  in  a  basic  task. 

There  are  three  main  factors  that  make  the  SIT  devices  less  attractive  in  RF  mobile 
applications.  First,  it  is  on  normally-on  characteristics,  second,  the  large  on-state  voltage 
drop,  and  third,  the  need  of  using  relatively  high  gate  bias  voltage  [1],  To  overcome  this 
problem,  neural  network  can  be  used.  We  obtain  much  SiC  SIT  data  of  properties  and  can 
be  used  for  training  neural  network  to  it. 

In  neural  network,  adequate  modeling  of  nonlinear  RF  devices  still  remains  a 
complex  but  an  avoidable  step  toward  the  achievement  of  a  good  design.  And  simulation 
boundary  is  very  flexible. 
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APPENDIX 
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A.  MEDICI  Program m  ing  Source  Examples. 


ASSIGN  NAME=vgtnin  C.VAL=0.0 

ASSIGN  NAME=vgstep  C.VAL=-2 

ASSIGN  NAME=nvg  C.VAL=1 1 

ASSIGN  NAME=vdmnx  C.VAL=200 

ASSIGN  NAME=Lsg  N.VALUE=0.5 

ASSIGN  NAME=RDTH  N.VALUE=50 

ASSIGN  NAME=MYSAT  N.VALUE=2.2E7 

ASSIGN  NAME=RS  N.VALUE=1 
ASSIGN  NAME=RD  N.VALUE=1 


LOOP  STEPS=6 

SASSIGN  NAME=MYSATN.VALUE=(L5E7  2,0e7  2.2c7  2.5c7  ) 
ASSIGN  NAME=TMPl  N.VALUE=(800  700  600  500  400  300) 


LOOP  STEPS=5 

ASSIGN  NAME=Lgatc  N.VALUE={0.05  0. 1  0.2  0.3  0.5) 
$  ASSIGN  NAME=Lgatc  N.VALUE=0.2 


LOOP  STEPS=5 

ASSIGN  NAME=RS  N.VALUE=(0.5  L  5  1050) 

\  LOOP  STEPS=5 

ASSIGN  NAME=RD  N.VALUE=(0.5  1  5  1050) 


LOOP  STEPS=5 
ASSIGN 
ASSIGN 
ASSIGN 
ASSIGN 
ASSIGN 


NAME=d 

NAME=d 

NAME=T 

NAME=t 

NAME=pilch 


N.VALUEKL75  1.5  L2  1,0  0.75) 
N.VALUE=L2 
N.VALUE=  1.5 
N_VALUE=  @d 
N.VALUE=  @d+@t 


CALL  FILE- 'acNd3_Thc_clcc_RSRD_ l.md" 


L.END 

L.END 

LEND 

L.END 

LEND 

STOP 
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«acNd3_Thc_c!cc_RSRD_  1 » 


TITLE  SIC  SIT  Temperature  Simulation 
COMMENT  finds  the  DC  and  small  signal  AC  parameters 


CALL  FlLE="Mobn.md" 


ASSIGN 

ASSIGN 

ASSIGN 


NAME=MYMAX 

NAME=MYMIN 

NAME=MYCON 


N.VALUE=@Mumax 

N.VALUE=@Mumin 

N.VALUE=<§NNrcf 


MESH  OUT.FILE="../.y../Mcsh/brcak_M@Str".mcsh" 


X.MESH 

X. MESH 

Y. MESH 
Y.MESH 
Y.MESH 
Y.MESH 


X.MIN  =0.0 

X. MIN  =  @d/2 

Y. MIN=-0.5 
DEPTH=(^Lsg 
DEPTH=@Lgatc 
DEPTH=<^Ldl 


WIDTH=@d/2  H  1=0.1 

WIDTH=@t/2 
DEPTH=0.5  HI  =0.1 

Hl=0.1  H2  =0.025 
HI  =0.025 
HI  =0.1  H2=0.25 


H  1=0.1 
H2=0.05 


REGION  NAME=SiC  SIC 

REGION  NAME=n_plus  SiC  SIC  Y.MAX=0.0 

REGION  NAME=air  INSULATO 

+  X.MIN=@d/2  X.MAX=@pitch/2  Y.MIN=-0.5  Y.MAX=@Lsg+@Lgalc/2 


ELECTR  NAME=sourcc  X.MAX=@d/2  TOP 

ELECTR  NAME=gate  VOID  X.MIN=@d/2  X.MAX=@pitch/2 

+  Y.MIN=@Lsg  Y.MAX=@Lsg+@Lgate 

ELECTR  NAME=drain  BOTTOM 

ELECTR  NAME=drainth  BOTTOM  THERMAL 

PROFILE  N-TYPE  UNIFORM  N.PEAK=@ND  OUT.  FlLE=^/../.. /Mesh/break JVgSt^pro” 

PROFILE  N-TYPE  UNIFORM  N.PEAK=1.2EI9 

+  Y.MAX=0.0 

CONTACT  NAME=gate  WORKFUNC=5.2  BARRIERL 

CONTACT  NAME=source  CON.RESI=@RS*(@d*lE-5) 

CONTACT  NAME=drain  CON.RES  I=((§jpitch/2)*(@RD*  1 E-5) 

CONTACT  NAME=drainth  R.THERM  A=(5)RDTH*  1 E3 

CALL  FILE="4h-SiC_pur_my.md" 

MATERIAL  INSULATO  PERMITTI=1.0 


MODEL  FERMID1R  INCOMPLE  ANALYTIC  FLDMOB  TEMPERAT=300 
SYMBOL  CARRIERS=1 
METHOD  CONT.STK 
SOLVE  INI 

REGRID  POTEN  IGNORE=air  RATIO=.l  MAX=I  SMOOTH=l 
+  IN.FlLE="../.y../Mcsh/brcak_"@Str".pro" 

+  OUT.  FILE-' .  J. ./.  ./McslVbrcak_M@Str" .  m  esh" 

REGRID  DOPING  LOG  IGNORE=air  RATIO=2 
+  IN.FILE="../.y../Mcsh/brcak_M@Str".pro" 

+  OUT.FILE=".y. ./ .  yMcsh/brcak_"(3!Str"  .mesh" 
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PL0T.2D  GRID  TITLE="Lg='Va!Lgate",d="@dM,Ldl=’,@Ldr.Nd="@ND 
+  BOUNDARY  DEPLETIO  FILL  SCALE  T.StZ£=0. 15 

$  +  DEVICECP  PLOT. OUT- 'psfiics/’VSjde  vice" .ps" 

PLOT.  1 D  DOPING  TITLE="Doping  at  x=0  "@dcvicc"  " 

+  X.START=0.0  X.END-0.0  Y.START-0.5  Y. END=@Lsg+@Lgale+®Ld  1 

+  LOGARITH  T.SIZE-0.2 

COMMENT  Compute  equilibrium  solution  as  initial  guess 
LOG  OUT .  F I  LE=" . /tc  mp/"@dc  vice",  log” 

SYMBOL  CARRIERS*  I  NEWTON 

METHOD  STACKS =20  CONT.STK 
SOLVE  V{gatc)=0.0  V(sourcc)=0,0  V(drain>=0.0 

REGRID  POTEN  lGNOR£=air  RATIO-05  MAX-1  SMOOTH*  l 


+  IN .  F  ILE= ./Mesh/break_"@S  tr" .  pro” 

+  OUT .  FI  LE  .  7.  A  ./Mesh/breaL_"@Str" .  m  esh" 


SYMBOL  CARRIERS- 1  NEWTON  ELECTRONS  COUP.LAT 

METHOD  STACKS-20  CONT.STK 

SOLVE  V(gate)=0.0  V(source)=0.0  V(drain)=0.0 

LOG  OUT .  FILE= /I  og/"(§idcvicc"  .log" 

SOLVE  V (gate)=@vgrmn  ELECTROD-gale  VSTEP=@vgstep  NSTEP=@nvg-l 

+  OUT.FILE=".7../../Gatc_Sol/"^Stj"_(K)"  SAVE.BIA 

CALL  FILE="biaspomi_T cmp_AC_anal . md " 
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«  Mobn.md» 

COMMENT  Mobility  Caculation  —  Depend  on  Temperature. 

ASSIGN  NAME=Mumax  N.VALUE=950*(@TMPl/300)**(-2.4) 

ASSIGN  NAME=Mumin  N.VALUE=40*(@TMP1/300)**<-0.5) 

ASSIGN  NAME=NNrcf  N.VALUE=2E17*(@TMP  1/300) 

ASSIGN  NAME=alphn  K  VALUE-0.76 
ASSIGN  NAME=Tbeta  N.VALUE=0.8 16 
ASSIGN  NAME=an  N.VALUE-4.27E-2 
ASSIGN  NAME=bb  N.  VALUE-98. 4 
ASSIGN  NAME=tt  N.VALUE=327 

ASSIGN  NAME=beta  N. VALUE=@Tbeta-H@afl*expf(@TMPl-@tt)/@bb) 


«  biaspoint_Tcmp_AC_onal-md» 

TITLE  Bias  point 

ASSIGN  NAME=vgnc  C.VAL--0.0 
ASSIGN  NAME=vdac  C.VAL=10 
ASSIGN  NAME=startso!  C.VAL=00 
C ALL  F IL  E="b  i  asp  lot_T cmp_AC  _anal_RSRD .  md" 

ASSIGN  NAME=vgac  C.VAL-0.0 
ASSIGN  NAME=vdac  C.VAL=25 
ASSIGN  NAME=startsol  CVAL=00 
CALL  FlLE="biasplot_Temp_AC_anal_RSRD.md" 

ASSIGN  NAME=vgac  C.VAL=-4.0 
ASSIGN  NAME=vdac  C.VAL-50 
ASSIGN  NAME=startsol  C.VAL-02 
C  ALL  F I  LE="  bi  aspl  ot_T  emp_AC  _an  a!_RS  RD.  md" 

ASSIGN  NAME=vgac  C.VAL=-6,0 
ASSIGN  NAME=vdac  C.VAL-75 
AS  S IGN  NAME=startsol  C.  VA  L=03 
C  ALL  F  ILE="  bi  aspl  ol_T empA  C_ona  l_R  S  RD.  m  d" 

ASSIGN  NAME=vgac  C.VAL—10.0 
ASSIGN  NAME=vdac  C.VAL=75 
ASSIGN  NAME=startsol  C.VAL=05 
CALL  FILE="biaspIot_Tcsip_AC_anal_RSRD.  md" 

ASSIGN  NAME=vgac  C.VAL=-4.0 
ASSIGN  NAME=vdac  C.VAL=35 
ASSIGN  NAME=startsol  C.VAL=03 
CALL  FILE="biaspiot_Tcmp_AC_anal_RSRD.mdn 
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«  biaspIot_Tcmp_AC_anal_RSRD.md» 

COMMENT  Small  signal  parameters  for  Vd=vdac ,  Vg=vgac 

TITLE  "VD~H@vdac"V,  VG-’^vgac'V  Str="@Strw  Lg-^Lgate"  d="@d”  Temp="@TMP  1  T.SIZE=0.2 

LOAD  IN.FILE=H.7../.7Gate_Sor@StrH_M@startsol 

METHOD  STACK=20  CONT.STK 

LOG  OUT.FILE«^./../.ytcmp/,,@dcvicc^log_NC,, 

SOLVE  V(sourcc)=0  V(gatc)=@vgac  V(drain)=@vdac 

PLOT. 2D  TITLE="VD="@vdac"V,  VG-"@vgacMv  Str="@Str"  Lg="@Lgatc"  d=M@d”Tcmp="@TMPr" 
T.SIZE-0.2 

+  BOUNDARY  DEPLETIO  FILL 
$+  DEVICE=0 

PLOT.OUT=".7.y.7psfiles/Temp  _"@StrM_,,@TMPrt_M@MYSAT'f_"@Lgate,,_"@d"_"@RS”J,@RD’’.ps” 

SOLVE  V(sourcc)=0.0  V(drain)=@vdac  V(gate)=@vgac 
+  AC.ANALY  FREQUENC-0.5E9  HI.FREQ  FSTEP-0.5E9  NFSTEP=20 

PLOT.  ID  IN. FILE*”. 7.7. 7tcmp/,,@deviccw. log_NCM  X.AXIS=freq 

+  Y.  AXIS ="G(gatc, gate)"  X.LOG  Y.LOG  TITLE="VD="@vdac"[VL  VG="@vgac"fv]  Str=M@Str"  Lg="@Lgatc" 
d="@d"  Tcmp=M@TMP  1 ""  T.SIZE=0.2 

40UT.FILE=M.7.7../Rcsults_acJ,@Str7TempJ,@TMPl"J,(2^YSAT"J’@Lgate"_"@d"J*@startsorj\2)vdac".ggg" 
S  PLOT.OUT=".7.7../psfiles/Tcmp_,,@Str"_,’@TMP  1  "_"@MYSAT"  J^Lgale'^^^NoContalor.ps" 

PLOT. ID  IN.FILE=M.7.7.7tcmp/t,@dcvicc".log_NC"  X.AXIS=frcq 

+  Y.AXIS=X}(drain,drainr  X.LOG  Y.LOG  TTTLE=MVD=,,@vdacKV,  VG="@vgac"v  Str="@Str"  Lg="@Lgatc" 
d="@d"Temp=3"@TMP  1  **  T.SIZE=0.2 
+ 

OUT.FILE=,,.7.7.7Rcsults_ac_',@Str,,/Temp_M@TMP^_n@MYSATt,_,,@LgateM_,,@d,'__"@sta^tso^’J,@vdac”.gddf, 
PLOT.  ID  IN.FILE-".  7.7.  ./temp/"@dcvice".Iog_NCM  >CAXIS=freq 

+  Y.AXIS="G<draw,gatc)M  X. LOG  Y.LOG  TITLE-” VD=w@vdacNV,  VG="@vgac"v  Str="@Str"  Lg="@Lgate” 
d**M@dnTemp=i,,@TMP  1 ""  T.SIZE=0.2 
+ 

OUT.FILE=".7.7.7Results_ac_,’@St^M/Temp_"@TMPI"_”@MYSAT,,_"@Lgatcn_"@d,'J’@sta^tso^,_,,@vdaC".gdg,, 
PLOT.  ID  IN.FILE=,,.7.7.7tcmpr@dc\icc’,.log-NC"  X.AXIS=fircq 

+  Y.AXIS="G(gate, drain)"  X.LOG  Y.LOG  TITLE="VD=,,@vdacMV,  VG-"@vgac"v  Str="@Str"  Lg="@Lgatc" 
d=,"@dwTempas"@TMP  1 ""  T.SIZE=0.2 
+ 

OUT.FILE«”.7../.7RcsuIts_ac_H@Str,,/Tcmp_"@TMPlM_,,@MYSAr_M@LgateV@d,,_r,@startsor,_,,@vdaCH.ggd" 
PLOT.  ID  IN . F ILE=" . 7. 7. 7tcmpr@dc  vice".  log^_NC "  X.AXIS-frcq 

+  Y.AXIS="G(drain,sourcc)"  X.LOG  Y.LOG  TITLE="VD="@vdac"V,  VG="@vgac"v  Str="@Str"  Lg="@Lgate" 
d=w@d"Temp=M@TMPrH  T.SIZE=0.2 
+ 

OUT.FILE*M.7.7.7Rcsults_ac__H@St^"/Temp_r,@TMPlV@MYSATV@Lgatc,,_,,@dV@sta^tsolM_"@vdac".gds,, 

PLOT.  ID  IN.FILE=".7../../tcmpr@devicc".log_NC"  X.AXIS=freq 

+  Y.AXIS="G(sourcc,drain)"  X.LOG  Y.LOG 

+  TITLE="VD="@vdac”V,  VG-"@vgac"v  Str»"@Str"  Lg="@Lgatc"  d="@d"Tcmp-"@TMPl"n  T.SIZE-0.2 
^UT.FILE*".y.y.7R(3idts_acJf@Str"/TempJ,@TMPl"J,@MYSAT"J,@Lgatc"J,@d"_,'@startsor,_"@vdac"  gsd" 

PLOT.  ID  IN.FILE=" .  7. ./.  7tcmp/"@dc  vice".log_NC  "  X.AXIS=freq 
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+  Y.  AXIS- 'C(gate,drain)"  X.LOG 

+  TITLE=’,VD='T@vdac,'V,  VG="@vgac"v  Str="@Str"  Lg="@Lgatc”  d=w@dwTcmp*H@TM PP"  T.SIZE=0.2 
40UT.FTLE="../. ./.  ./Rcsults_ac Jf  @StrH/Tcmp_M@TMP  1  M_M@MY  S  AT"  J,@Lgatc"_,’@d0_M@startsol,,_f,@vdacH.cgd" 

PLOT.  1 D  IN.FILE=M../J../tcmpr@devicc,Mog_NCn  XLAXIS-frcq 

4-  Y.AXIS="C(gate,  source)”  X.LOG 

+  TITLE=" VD="@vdacMV,  VG=M@vgac"v  Slr=H@Sl r"  Lg="@LgatcM  d«M@d,Tcmp=M@TMPlm,  T.SLZE=0.2 
+  OUT.FILE=M.7../../Results_acJ,@StrM/Tcmp_”@TMPr_”@MYSAT"J,@Lgatcn__M@d,,_M@startsor,Jt@vdac,,.cgsM 

PLOT.  1 D  IN.FILE=\./.7.7tOTip/"@devic^logJ^C"  X.  AXIS=freq 

+  Y.AXIS=^t'C(drain,sou^ce),,  X.LOG 

4-  TITLE=”  VD="@vdac"V,  VG="@vgac*'v  Str=”@Str"  Lg=”@Lgate"  d="@d"Temp=M@TMP  1 ""  T.SIZE=0.2 
+  OUT.FILE=l'^./../../Rcsulls_ac_"@Str,'/Temp_,,@TMP^_M@MYSAT"J'@Lgate,,_,,@d'’_,t@sta^tso^'Jf@vdac,,.cds,, 

PLOT.  ID  IN.FILE=H../../../tcmpr@deYicc,,.log^NCtt  X.AXIS=freq 

+  Y.AXIS="C(gate,gate)M  X.LOG 

+  TITLED  VD*"@vdac"V,  VG="@vgac"v  Str*”@Str*  Lg=f,@Lgate"  d=,,@dT,Tcmp=M@TMPr,M  T.SIZE=0.2 

4- 

OUT.FILE=’'.  7.7.  ./Results_ac_t,@Str"/Tcmp_,’@TMP  l,,_"@MYSAT^_,'@Lgalc’,_M@d,,_,,@sla^tso^’_,,@vdac".cgg,, 

PLOT.  ID  IN.FILEs=M.7.7.7lemp/t,@dcvice".Iog_NCM  X.  AXIS-freq 

+  Y.AXIS=,,C(soiircc,drain)f'  X.LOG 

4-  TlTLE="VD-"@vdac,'V,  VG=”@vgac"v  Str="@Str"  Lg="@Lgate"  d=,,@d,,Tcmi>=w@TMPrM  T.SIZE=0.2 
4-OUT.FILE=n../../.7RcsuUs_ac_,,@Str,,/Temp_"@TMPr_"@MYSAT,,_"@Lgate,,>_"@d',_,,@startsol't_M@vdac".csdn 

PLOT.  ID  IN.FILE=^7.7.7tcmp/,,@deviceMogL_NCl,  X.AXIS-frcq 

+  Y.AXIS=HC(gatc,sourcc)w  X.LOG  Y.LOG  ABSOLUTE 

4<>UT.  FILE=,,../.7../ResuIts_acJ,@Str,,/Temp_M@TMP  lM_M@MYSAT,,J'@Lgate,,_,,@d,,_M@startsor’_H@vdacn.cgs  1 

n 

PLOT.  ID  IN.FILE=’,../.7.7temp/n@dcvice',.!ogL_NC',  X.AXIS=freq 

+  Y.  AXIS=" Y(gatc,drain)M  X.LOG  Y.LOG  ABSOLUTE 

4-  Tm.E="VD="@vdacMV,  VG=,,@vgac,,v  Str=H@Str"  Lg="@LgateM  d=ft@dMTcmp=w@TMPr,w  T.SIZE=0.2 

4- 

OUT.FILE=".7.y.iResults_ac_"@Str"/Temp_"@TMPr_"@MYSAT"_”@Lgate"_"@d"_"@staitsol"_H@vdac".yl2" 

PLOT.  ID  IN.FILE="./y../terapr@device".log_NC"  X.AXIS=frcq 

+  Y.AXIS="Y(drain,drain)"  X.LOG  Y.LOG  ABSOLUTE 

+  TITLE- 'VD="@vdac"V,  VG="@vgac"v  Str="@Str"  Lg="@Lgatc”  d="@d"Temp="@TMP  1""  T.SIZE=0.2 
+ 

OUT.FILE="./y../Rcsults_ac_H@Str"/Teinp_"@TMPr_"@MYSAT"_"@Lgate"_"@d"_"@startsor_”@vdac".y22" 

PLOT.1D  IN.FlLE=”.y../.ytcmpr@devicc,,.logL.NC"  X.AXIS=freq 

+  Y.  AX1S="  Y (gate.gate)"  X.LOG  Y.LOG  ABSOLUTE 

+  TITLE- 'VD=”@vdac"V,  VG="@vgacHv  Str="@Str"  Lg="@Lgate"  d=B@d"Temp=”@TMPI""  T.SIZE=0.2 

4 

OUT.FILE=".y.y../Rcsults_ac_"@Str"/Tcmp_"@TMP  l”_"@MYSAT"_,’@Lgatc,,_"@d"_"@startsol"_"@vdac''.y  1 1  ” 
PLOT.  1 D  IN .  FILE=" . ./.  J.  ,/tempr@de  vice"  ,log_NC"  X.AXlS=freq 

+  Y.AXIS="Y(drain>gatc)"  X.LOG  Y.LOG  ABSOLUTE 

+  nTLE="VD="@vdac"V,  VG="@vgac”v  Sto="@Slr"  Lg=”@Lgate"  d="@d"Temp=H@TMPr"  T.SIZE=0.2 

+ 

OUT.FILE- \7../.  7Results_ac_M@Str7Temp_w@TMP  1  ,’-,,@MYS  AT”  J,@Lgate<'_*,@dM_M@startsol"__M@vdac,’.y2 1 " 

PLOT.  ID  IN . F ILE=" . 7.7. 7tempr@device" .log_NC  "  X.  AXIS=frcq 

4-  Y.AXIS-"  Y(sourcc,dnuD)w  X.LOG  Y.LOG  ABSOLUTE 

4-  TITLE==’'VD="@vdac"V,  VG=tt@vgac"v  Str="@Str"  Lg=,r@Lgate"  d="@ dTemp="@TMPr"  T.SIZE-0.2 
4* 

OUT.FILE88".  7.7.  ./Results__acJ,@Str"/Temp_"@TMP  YSAT,_,,@Lgate,'_"@d,,_M@sta^tsol,,_',@vdac".ysd,, 

STOP 
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B.  MATLAB  Example  Sources 


«Tcmp_frcq_Str_alL_calrm» 
close  all 
clear  all 
format  long  g 

%  Structure  type 
%for  Sir*  1:4 
Str^l; 

LLg  =  zeros(6,4i4); 

DD1  *  zeros(6,4,4); 

CCgs  =  zeros(6,4,4); 

C Cgd  -  zeros(6,4T4); 

FFt  **  zcros{6J4,4)i 

%  Concentration  Nd 
o^Slr; 

%  Velocity  drift  saturation 

p-5; 

rs=2;rd^2; 

%m=l; 

%for  rs=l  :4 

if  Sir^l 

d- [1.75  1.5  L2  1.0]; 
elseif  5tr=2 
d=[L2  10.8  0.6]; 
clscif  $lr^=3 

[0.8  0.7  0,6  0.5  0.4]; 

else 

d-  [0.6  0.5  0.4  0.3  0.25]; 
end 

V  5ai=f  ’  1  <  0E+071; '  1 15  E+0772 ,0E+0T;  '2 .2E+OT;7,5E+0  7']; 

L*  [*0,50000';  ’0.20000';  '0. 10000';  '5.0c-021; 

Tempr-fSOO'^OGYSOO'iW^GOYSOOl; 

Nd-f5.0e+I5*;  'L0&H61;  ’S-Oc+ieVZOc+ir]; 

RS=[0.5  1  5  10]; 

RD=[0,5  !  5  101; 

fortt^lrfi 

i=0; 

forn**l:5 
for  m=I:5 
i-i+1; 

ifn=5 

FH[‘ .7..fTcmpr/Results_ac_’  nuni2slr(Sfr)  YTcmpJ  Tempr(tt,:)  rJ  Vsat(p,:)  *J  L(n,;)  'J  num2str(d(m)) 
'_NoContator_f  num2sti(RS(rs))  \y  1 1*]); 

F2=<[1../../Tcmpr/Rjcsults_ac_'  nuni2str(Str)  '/Temp  J  Tcmpr(tt,:)  *J  Vsat(pT:)  *J  L(n,:)  *J  nmn2str(d(m)) 
4_NoContator'  mim2str<RS(rs))  l.yl2*J); 

F3^f  ./.yTempr/R£sults_ac_f  mim2str(Str)  VTcmpJ  Tcmprftt,:) Vsat(p,:)  rJ  L(n,:)  *J  num2slr(d(m)) 
*_NoContatorJ  num2str(RS(rs))  '_y21T); 

F4=<[l  /.frcmpr/Results_ac_l  mun2str(Str)  VTcinpJ  Tcmp^tt,:)  *J  Vsat(p,:)  'J  L(n,:)  *J  num2slr(d(m» 
’NoContatorJ  num2sti(RS(ra))  Vy22']); 
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F5Kr-././rempr/Results_acJ  num2str(Slr)  '/TempJ  Tempr(tt,:) Vsat(p,:)  'J  L(n,:) 'J  num2str(d(m)) 
JNoContatorJ  num2stifRS(rs))  '.cgd']); 

F6=([\./.  /Tempr/Resultsac J  num2str(Str)  '/TempJ  Tempr(tt,:)  'J  Vsatfp,:) L(n,:) 'J  num2str(d(ni)) 
JNoContatorJ  num2str(RS(rs))  '.cgs']); 

F7=([,../../Tcmpr/Resalls_ac_’  num2str(Str)  '/TempJ  Tempr(tt,:) J'  Vsat(p,:)  J*  L(n,:)  'J  num2str(d(ra)) 
JNoContatorJ  num2str(RS(rs))  '.gds']); 
else 

Fl^p./VTcmpr/ResultsacJ  num2str(Str)  '/TempJ  Tcmpr(tt,:)  'J  Vsat(p,:)  *J  num2str(str2num(L(ii, : )))  J 
num2str(d(m))  JNoContatorJ  num2str(RS(rs))  '.yl  I*]); 

F2={[,.y../Tempr/Results_acJ  num2str(Str)  '/TempJ  Tcmprftt,:)  'J  Vsat(p,:)  'J  num2str(str2num(L(n7:)))  'J 
num2str(d(m))  JNoContatorJ  num2str(RS(rs))  '.yl2']); 

F3=(['../../Tempr/ResuIts_acJ  num2str(Str)  '/TempJ  Tempr(tt,:)  'J  Vsat(p,:>  'J  num2str(str2num(L(n,:)))  J' 
num2str(d(m))  JNoContatorJ  num2str(RS(rs))  '.y21r]); 

F4={[,../.7Tcmpr/Results_acJ  num2str(Str)  ’/TempJ  Temprftt,:)  J'  Vsat(p,;)  J  num2str(str2nuin(L(n,:)))  'J 
num2str(d(m))  JNoContatorJ  num2str(RS(rs))  \y22*]); 

F5=<[,./7Tempr/RcsuIts_acJ  num2str(Str)  '/TempJ  Temprftt,:)  'J  Vsat(p,:)  'J  num2str(str2niim(L(n,:)))  J' 
num2str(d(m))  JNoContatorJ  num2str(RS(rs))  \cgd']); 

F6*((’../../Tcmpr/Rcsiilts_acJ  num2str(Str)  ’/TempJ  Temprftt,:)  'J  Vsat(p,:)  'J  num2str(str2num(L(n,:)))  'J 
num2str(d(m))  JNoContatorJ  num2str(RS(rs))  '.cgs']); 

F7=(r../../Tempr/ResuIts_acJ  num2str(Str)  '/TempJ  Temprftt,:)  'J  Vsat(p,:)  'J  num2str(str2num(L(n,:)))  'J 
num2str(d(m))  JNoContatorJ  num2str(RS(rs))  '.gds']); 
end 

[fid  1  ,err  1  ]=fopen(F  1  f'r'); 

[fid2,crr  1  ]=fopen(F2,'r’); 

[fid3,errl]=fopcn(F3,V); 

[fid4,eir  l]=fopen(F4,'r'); 

[fid5,errl]*fopen(F5,'r'); 

[fid6,crr  1  ]  =fopen(F  6,'r1); 

[fid7terr  1  j=fopen(F7,'r’); 

Sll-fccanf(fidl/%s’,15); 

5 12  =  fscanfl[fid2,'%s',l5); 

513  =  fscanF[fid3,'%s',15); 

514  =  fscanf(fid4,'%s',15); 

5 1 5  -  fscanfffidS/Vos',  1 5); 

516  =  6canfl[fid6I>%s,,15); 

S17«fscanfi[fid7,,%8',15); 

541  -  fscanfffidl,  ’%16e',[2  inf]); 

542  ~  fscanf(fid2,  '%16e',[2  inf]); 

543  «  fscanf(fid3,  '%\6c\[2  inf]); 

544  =  fccanfl(fid4,  K,/«16e\t2  inf]); 

545  =  fscanfl[fid5,  *%16e',[2  inf]); 

546  =  fscanf(fid6,  *%16e\[2  inf]); 

547  =  fscanf(fid7,  '%16c',[2  inf]); 

fclosc(fidl);fclose(fid2);fclose(fid3);fclose(fid4);fclose(fid5); 
fdose(fid6)  ;fclose(fid7) ; 

S21=('Frcquency  GHz'); 

S3  l=('Condiictancc(pF/inm)'); 

yl  =  S41(2,:); 
y2  =  S42(2,:); 
y3  =  S43{2,:); 
y4  =  S44(2,:); 
y5  =  S45(2,:); 
y6  =  S46(2,:); 
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y7  «  S47(2,:); 

xl  =  (S41(l ,:)) 
x2  •  (S42(l,:)) 
x3  =  (S43(l ,:)) 
x4  =  (S44(l ,:)) 
x5«(S45(l,:)) 
x6- (846(1,:)) 
x7  =  (S47(1,:)) 

f~xl; 

w=f.*2*pi; 

yi=yl-y2; 

yo=y3-y2; 

yr=-y2; 

yf=y3-y2; 


Cgd=yr./\v; 

Cgs=-yi./w; 

Y 1 2=complex(0,-w.  *y5); 

gml=y3+w.*y5; 

gm2=abs(y0; 

gds=abs(yo); 

n5=length(y5); 

Lg=str2num(L(n, :  ));Lg; 

DD=num2str(d(m));DD; 

nn-length(gml); 

gm  1 1  *=(sum(gm  l)/nn)*  1  c6; 

gm22=(sura(gm2)/im)  *  1  e6; 

n5=length(y5); 

Cgdl  1  =(sum(abs(y5(2:n5)))/(nn- 1 )).  *  l  e  1 5; 

Cgd22  =(sum(abs(Cgd(2  :n5)))/(nn- 1 )). *  le  1 5 ; 

Cgsl  1  =(sum(abs(y6(2  :n5)))/(nn- 1 )).  *  1  c  1 5 ; 

Cgs22  Ksum(abs(Cgs(2:n5)))/(nn- 1  )).*  lcl  5; 

Ftl=gml  l/(2*pi*(Cgsl  i+Cgdl  1)); 

Ft2=gm22/(2*pi^Cgs22+Cgd22)); 

LLgl(tUvm)  =  Lg; 

DDI  l(tMMn)  *  str2num(DD); 

GGml(ll,iyn)  =  gm22; 

CCgsl(tt,n,m)  -  Cgs22; 

CCgd  K^m)  =  Cgd22; 

FFtl(tMMn)  =  Ftl; 

VVsatftMyn)  *  slr2num(Vsat(p,:)); 

%dispC . . . .  ■  *); 

end 

end 

end 

fclose  all; 

Filename  =  (rdi/Research/DaTa/Frequency/NoContact/fre^StrJ  num2str(Str)  '.mat']); 
save(Filename,  'Wsat',  'LLgl',  ’DDll’,  'GGml’,  ’CCgsl’,  ’CCgdl\  ’FFtlVRS*); 

%end 

%end 

diSp(t— . . . -  — -  ==«»»); 

disp(’===^=^^— Strjill_Cal==*^--^ — =T); 
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«  Tcmp_freq_Str.m» 


% 

%  Include  contact  Resistor 
% 


close  all 
clear  all 
format  long  g 


%  d  is  3  steps 
Lg=zcros(l,5); 

Tempr=r300,;,400,;,500,;,600,;700,;,800’]; 

Vsat=['1.0E+07,;,1.5E+07,;,2.0E-K)7f;,2.2E+07,;,2.5E407t]* 

RS=[0.5  1  5  10  50]; 

REK0.5  1  5  10  50]; 

Nd=[’5.0e+15';  U.Oe-t-l^;  ,5.0e+16';,2.0e+17T]; 

Bias~[.0, 10;  0.0, 25;  -4, 50;  -6.0,  75;  -10.0,  75;  -4.0,  35]; 

vgac=Bias(:,l)'; 

vdac=Bias(:r2)'; 

Sti=4; 

BA=1; 

rs“2;rd=2; 

if  Str=l 

d=  [1.75  1.5  1.2  1.0  0.75]; 
elseifStr=2 

d=[1.4  L2  10.8  0.6]; 
elseif  Str==3 

d=  [0.8  0.7  0.6  0.5  0.4]; 
else 

d=  [0.6  0.5  0.4  0.3  0.25]; 
end 

fi-i; 

ii=l; 

fortt=l:6 

Filename  =  (['g:\Temprt\Data\lreq_Str_'  num2str(Str) '  '  num2str(Tempr(tt,:))  ’  '  num2sti(BA)  V  num2str(RS(rs)) 
num2slr(RD(rd))  '.mat']); 
load(Filcnamc,  ’-mat'); 

Lg  =LLg(:,l)’; 
d  -DD1(1,:); 

evalc([7t_’  num2str(Tenipr(tt,:))  -FFt']); 
evalcCj’Gml  J  num2str(Tempr{tt,:)) '  =GGm’]); 
evalcff'CgslJ  num2str(Tempr(tt,:)) '  =CCgs']); 
evalc([’Cgdl  J  num2str(Tenipr(tt,:)) '  =CCgdl); 
evalc([’Vsat  =Wsat’J); 


%RS  =  RS.*10*-5; 
N=40; 
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Nl=40; 

nx  1  =(d(  1 )  -  d(5))/N; 

nyl=(Lg(l).Lg(5))/N; 

ny2  =  (RS(1,5)-RS(1,1))/N1; 


[xl,yl]  *  meshgrid(  d(l):-nxl:d(5),  Lg(l):-nyl:Lg(5)); 

[x2,y2] -  mcshgrid(  Lg(l):-nyl:Lg(5),  RS(I,l):ny2:RS(l,5)); 

[x3,y3]  =  meshgrid(  d(l):-nxl  :dC5X  RS(I,l):ny2:RS(l,5)); 

evalctfzilj  num2str(Te!npr(tt,:)) '  *=  interp2(d,  Lg,  FtJ  num2str(Tenipr(tt,:)) xl,  yl)’  ]); 
%  change  data  structure 

%  cvalcflTtl  J  num2str(m)  JRS(:,rs)  »  FtJ  num2str(m)  JTJ  num2str(rs)  '(:,3)’]); 

%  eva!c([TFt2J  num2slr(m)  JRS(rs,:)  =  FtJ  num2str(m)  JTJ  num2str(rs)  '(3, :)*]); 


figurc(I); 

cvalcd^surfc^xl,  yl,  adlj  num2str(Tcmpr(tt,:))  ')']); 

%  evalcft'fXXl  1,  YY1 1]  =  max(/ir  num2str(n>)  J]); 

%  evalc([’ZZi=zir  num2str(rs)l); 

%  tcxt(XXll,YYl  1, ZZi,  ['Maxft-,  num2str(RS(rs))  ]); 
hold  on 

xlabel('Mesa  Width  [{\it\mum}Jf); 
ylabcl('Gatc  length  [{\it\mum}]'); 
zlabelC  Cut-off  Frequency  [{\itGHz}]'); 

colorbar 
axis  tight 
axis  vis3d 
shading  interp 

figure(300); 

subplot(2,3j5) 

evalc(rsiirfc(xl,  yl,  zil  J  num2str(Tcmpr(tt,:))  7]); 
hold  on 

if  fi==l 

xlabel('Mesa  Width  [{\it\imun}]'); 

ylabclCGatc  length  [{\it\mum}]'); 

zlabelC  Cut-off  Frequency  [{\itGHz})*); 

titled  Temperature  « '  num2str(Tempr(tt,:))  ^iK^l'Cclor','!4) 

end 

if  fi  =2 

title<[  num2str(Tempr(tt,:))  \[{\itK}]VColorJr') 
end 

if  fi  — ^3 

title([  num2str(Tempr(tt,:)) {\itK }  ]']  JColor^JO 
end 

if  fi  =4 

title([  num2str(Tcnipr(tt,:))  ,,[{\itK}],],,CoIor,,V) 
end 

if  fi  «=»5 

titlc([  num2str(T empr^tt, : ))  ’,[{\itK}  f], 'Color',  V) 
end 

colorbar 
axis  tight 
axis  vi$3d 
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shading  intcrp 

cvalc(['zi4'  num2str(Tempr(U,;)) 1  =  interp2(d,  Lg,  FtJ  num2str(Tcmpr(tt,:))  x I,  y  1)'  ]); 
evalc^ztf'  num2str(Tempr(tt,:)) '  =  interp2(d ,  Lg,  GmlJ  num2str(Tcmpr(tt,:)) xl,  yl)'  ]); 
evalc^'ztf’  num2str(Tempr(tt,:)) ' -  interp2(d,  Lg,  CgslJ  num2slr(Tempr(tt,:))  \  xl ,  yl)'  ]); 
evalc(['zi7'  num2str(Tcmpr(U,:)) '  =  interp2(d,  Lg,  CgdlJ  num2stj(Tempr<tt,:))  xl,  yl)'  ]); 

figurc(60); 

evalc^'surfcfal,  yl,  zi7'  num2str(Tempr(tt,:))  *)']); 
hold  on 

xlabcl(’Mcsa  Width  [\mum]'); 
ylabcl(’Gate  length  [\mum]’); 
zlabcl('  Gate-Drain  Capacitance  [pF/mm]'); 

%CitIe<[a  Temperature  = '  num2str(Tctnpi(tt,:))  ,[KJ'l,,Color,,,rr); 

colorbar 

axis  tight 

%axis([0.05  0.5  (0.5*10A-5)  (10*10A-5)  6  14]); 
axis  vis3d 
shading  intcrp 

figure(70); 

cvalcd'surfcCxl,  yl,  zi6'  num2str(Tempr(tt,:)) *)D; 
hold  on 

xlabel(’Mesa  Width  [\mum]'); 
ylabcl('Gate  length  [\mum]'); 
zlabelC  Gate-Source  Capacitance  [pF/mm]'); 

Votitlcd*  Temperature  =  '  num2str(Tcnipr(tt,:))  '[K]f], 'Color1,  V); 

colorbar 

axis  tight 

%axis([0,05  0.5  (0.5*  1 0M)  (10*10^5)  6  14]); 
axis  vis  3d 
shading  intcrp 

figurc(71); 

subp!ot(2,3,ii) 

evaicd'siirfc(xl,  yl,  zi6'  num2str(Tcmpr(tt,:)) ')']); 
hold  on 

xlabelCMesa  Width  [\rnum]'); 
ylabel(’Gate  length  [\mum]'); 
zlabcl(’  Gate-Source  Capacitance  [pF/mm]'); 

%title([f  Temperature  = 1  num2str(Tempr(tt,:))  '[K]'],’Color',Y); 

colorbar 

axis  tight 

%axis([0.05  0.5  (0.5*10A-5)  (10*10A-5)  6  14]); 
axis  vis3d 
shading  interp 


figure(75); 

evalc(['surfc(xl,  yl,  zi7'  num2str(Tcmpr(tt,:))  ’)']); 
hold  on 

xlabcl('Mesa  Width  [\mum]’); 
ylabelfGate  length  [\mum]'); 
zlabelC  Gate-Drain  Capacitance  [pF/mm]'); 

%titlc(['  Temperature  = '  num2str(Tempr(tt,:))  '[K]'],,Color,,Y); 

colorbar 

axis  tight 

%axis((0.05  0.5  (0.5*10A-5)  (10*10A-5)  6  14]); 
axis  vis3d 
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shading  interp 


figurc(76); 

subplot(2,3,ii) 

evatc^surfcCxl,  yl,  zi7'  num2str(Tempr(tt,:)) '/]); 
hold  on 

xlabelCMesa  Width  [\imun]'); 
ylabel('Gate  length  (\mum]r); 
zlabcl('  Gate-Drain  Capacitance  [pF/mm]'); 

%titlc({,  Temperature  = '  num2str(Tenipr(tt,:))  ,[K]f] , ’Color1, V); 

oolorbar 

axis  tight 

%axis([0.05  0.5  (0.5M0A-5)  (10*i0A-5)  6  14]); 
axis  vis3d 
shading  interp 

figurc(80); 

evalc([’surfc(xl,  yl,  zt5'  num2str(Tempr(tt,:))  *)’]); 
hold  on 

xlabelfMesa  Width  [Amum]1); 
ylabel('Gate  lengtli  [\mum]'); 
zlabel('  Transconductance  [mS/mm]'); 

%titlc(f  Temperature  “  ’  num2str(Tcmpr(tt,:))  '[Kn/Color’/r’); 

colorbar 

axis  tight 

axis  vis3d 

shading  interp 

figure<81); 

subplot(2,3,ii) 

evalc(['$urfc(xl,  yl,  zi5’  num2slr(Tempr(tt>:)) ')’]); 
if  ii=l 

xlabel('Me$a  Width  [\mum]'); 
ylabel(Xjate  length  [\mum]'), 
zlabclC  Transconductancc  [mS/mm]'); 

titled'  Temperature  = '  num2str(Tcmpi<a,:))  MKn/ColorVr') 

else 

titled  nun^strCTempKtl,:))  ,,[K]r],,Color,,,0 
end 

%Utled'  Temperature  = '  num2str(Tcmpr(tt,:))  '[K]'],’Color\V); 

colorbar 

axis  tight 

axis  vis3d 

shading  interp 

figurc(90-Hi); 

evalc(['surfc(xl,  yl,  zi4’  num2str(Tempr(tt>:))  01); 

xlabel('Mcsa  Width  [\mum]'); 

ylabelCGatc  length  [\mum]'); 

zlabel('  Cut-off  Frequency  [GHz]'); 

title<r  Temperature  = '  num2str(Tempr(tt,:))  ’[K]'], 'Color',  V); 

colorbar 

axis  tight 

%axis«0.05  0.5  (0.5*10A-5)  (10*i0A-5)  6  14]); 
axis  vis3d 
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shading  intcrp 

iHi+1; 

end 


fclose  all; 
hold  off; 
dispO™5 
dispC== 
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Simulation  Results  for  I-V  Characteristics 
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Figure  C.l  I-V  curves  for  different  values  of  Lg  and  2a,  with  Nd  -  5-101;>  [cm']  and  300  [K\. 
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Figure  C.2  I-V  curves  for  different  values  of Lg and  2a,  with ND  =  5- 10 15  [cm3]  and  400  [K\. 
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Figure  C.3  I-V  curves  for  different  values  ofZ,s  and  2a,  withJVD  =  S-IO15  [cm3]  and  500  [AT). 
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Figure  C.4  I-V  curves  for  different  values  of  Lg  and  2a,  with  Nd  =  5-  10l3  [cm'3]  and  600  [K\. 
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Figure  C.5  I-V  curves  for  different  values  of  Lg  and  2a,  with  ND  =  5-1015  [cm'3]  and  700  [K], 
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Figure  C.6  I-V  curves  for  different  values  of  Lg  and  2a,  with  ND  =  5- 1015  [cm'3]  and  800  [AT], 
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Figure  C.7 1-V  curves  for  different  values  ofL,  and  2a,  with  No  =  MO16  [cm3]  and  300  [AT|. 
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Figure  C.8  I-V  curves  for  different  values  ofZ,g  and  2a,  with  No  =  MO16  [cm3]  and  400  [AT], 
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Figure  C.9  I-V  curves  for  different  values  of  Lg  and  2a,  with  Nd  =  1- 1016  [cm3]  and  500  [K\. 
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Figure  C.10  I-V  curves  for  different  values  of  Lg  and  2a,  with  No  =  MO16  [cm3]  and  600  [A!]. 
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Figure  C.  1 1  I-V  curves  for  different  values  of  Lg  and  2a,  with  ND  =  1  - 1 016  [cm'3]  and  700  [AT]. 
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Figure  C.  12  1-V  curves  for  different  values  of  Lg  and  2a,  with  No  —  1  1016  [cm1]  and  800  [X]. 
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Figure  C.13  I-V  curves  for  different  values  of  Lg  and  2a,  with  No  =  5  •  1 0 ' 6  [cm3]  and  300  [AT], 
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Figure  C.  1 4 1-V  curves  for  different  values  of  Lg  and  2a,  with  Afo  =  5- 1016  [cm'J]  and  400  [A], 
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Figure  C.15  I-V  curves  for  different  values  of  Lg  and  2a,  with  ND  =  5-1016  [cm3]  and  500  [£]. 
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Figure  C.16 1-V  curves  for  different  values  of  Lg  and  2a,  with  Afo  =  5-1016  [cm'3]  and  600  [K\. 
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Figure  C.17 1-V  curves  for  different  values  of  Lg  and  2a,  with  Nd~  5-1 016  [cm3]  and  700  [AT]. 
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Figure  C.18 1-V  curves  for  different  values  of  Lg  and  2a,  with  ND  =  5-1016  [cm3]  and  800  [K\. 
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Figure  C.19  I-V  curves  for  different  values  of  Lg  and  2a,  with  No  =  2- 1 0 17  [cm'3]  and  300  [AT]. 
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Figure  C.20 1-V  curves  for  different  values  of  Lg  and  2a,  with  No  =  2- 1 017  [cm3]  and  400  [AT]. 
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Figure  C.21 1-V  curves  for  different  values  of  Lg  and  2a,  with  No  =  21017  [cm3]  and  500  [AT|. 


109 


100  150  200 


50  1  00  ISO 

Drain  Voltaga  [V] 


. 

2= 

100  150  200 


Drain  VolUga  [V] 


50  100  190 

Drain  VolUga  (V] 


Oraln  VolUga  [VI 


I 

a 
I 

50  100  ISO  200 

Oraln  VolUga  [V] 


60  100  ISO 

Drain  VolUga  [V] 


Ij 

| . 

L . j 

. 

. 1 

L 

i 

i  1 

: . ] . . 

L . :i 

W'\ 

_ J 

1  1 

1’  i 

SO  100  ISO  200 
Drain  Voltaga  [V] 


SO  100  ISO  200 
Maaa  Width  d-  0.0  (,m] 


BO  100  ISO  200 
Maaa  Width  d«0.S  Urn] 


SO  100  150  200 

Mata  Width  d-  0.4 


Figure  C.22 1-V  curves  for  different  values  of  Lg  and  2a,  with  No  =  2-1017  [cm'3]  and  600  [XT], 
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Figure  C. 23  I-V  curves  for  different  values  of  Lg  and  2a,  with  ND  =  2-1017  [cm's]  and  700  [K\. 
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Figure  C.24 1-V  curves  for  different  values  of  Lg  and  2a,  with  ND  =  2-1017  [cm3]  and  800  [A], 
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D.  Conversion  from  Y-parameter  to  S-parameter 
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E. 


Simulation  Results  for  Frequency  and  Transconductance  Characteristics 


Gate  length  [^m]  Mesa  Width  fcm] 


Figure  E.l  Cut-off  frequency  for  different  values  of  Lg  and  2a,  with  No  =  5*10,s  [cm 
J|,  from  top  to  bottom,  300,  400,  500,  600, 700, 800  [A],  respectively. 
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Figure  E.2  Cut-off  frequency  for  different  values  of  Le  and  2a,  with  jV0  =  1*1016  [cm 
3\,  from  top  to  bottom,  300, 400,  500, 600,  700,  800  [A],  respectively. 
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Figure  E.3  Cut-off  frequency  for  different  values  of  Lg  and  2a,  with  ISfD  =  5*1016  [cm 
*|,  from  top  to  bottom,  300, 400,  500,  600,  700,  800  [A],  respectively. 
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Figure  E.4  Cut-off  frequency  for  different  values  of  Lg  and  2a,  with  No  —  2*1017  [cm 
J),  from  top  to  bottom,  300, 400, 500, 600,  700, 800  [A],  respectively. 
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Figure  E.5  Transconductance  for  different  values  of  Lg  and  2a,  with  IWD  =  5*1015 

[cm3 1,  from  top  to  bottom,  300, 400,  500,  600,  700,  800  [AT],  respectively. 


Gate  length  [pm]  Mesa  Width  [pm] 

Figure  E.6  Transconductance  for  different  values  of  Lg  and  2a,  with  Np  =  l-lO16 

[cm3 1,  from  top  to  bottom,  300,  400, 500,  600,  700,  800  [A],  respectively. 
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Figure  E.7  Transconductance  for  different  values  of  Lg  and  2a,  with  No  =  5-101 

[cm3 1,  from  top  to  bottom,  300,  400,  500, 600,  700,  800  [Z],  respectively. 
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Figure  E.8  Transconductance  for  different  values  of  Lg  and  2a,  with  No  =  2-1 017 

[cm%  from  top  to  bottom,  300,  400,  500, 600,  700, 800  [Jf|,  respectively. 
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